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ABSTRACT 

Aims. Despite their importance to a number of astrophysical fields, the lifecycles of very massive stars are still poorly defined. In order 
, to address this shortcoming, we present a detailed quantitative study of the physical properties of four early-B hypergiants (BHGs) 

of spectral type Bl-4 Ia + ; Cyg OB2 #12, Sco, HD 190603 and BP Cru. These are combined with an analysis of their long-term 
spectroscopic and photometric behaviour in order to determine their evolutionary status. 

Methods. Quantitative analysis of UV-radio photometric and spectroscopic datasets was undertaken with a non-LTE model atmo- 
sphere code in order to derive physical parameters for comparison with apparently closely related objects, such as B supergiants 
(BSGs) and luminous blue variables (LBVs), and theoretical evolutionary predictions. 

Results. The long-term photospheric and spectroscopic datasets compiled for the early-B HGs revealed that they are remarkably 
1 stable over long periods (>40 yrs), with the possible exception of ( l Sco prior to the 20'* century; in contrast to the typical ex- 

cursions that characterise LBVs. Quantitative analysis of Sco, HD 190603 and BP Cru yielded physical properties intermediate 
between BSGs and LBVs; we therefore suggest that BHGs are the immediate descendants and progenitors (respectively) of such 
stars, for initial masses in the range ~ 30 - 60M G . Comparison of the properties of ^' Sco with the stellar population of its host clus- 
ter/association NGC 6231/Sco OBI provides further support for such an evolutionary scenario. In contrast, while the wind properties 
of Cyg OB2 #12 are consistent with this hypothesis, the combination of extreme luminosity and spectroscopic mass (~ 110M Q ) and 
comparatively low temperature means it cannot be accommodated in such a scheme. Likewise, despite its co-location with several 
LBVs above the Humphreys-Davidson (HD) limit, the lack of long term variability and its unevolved chemistry apparently excludes 
such an identification. Since such massive stars are not expected to evolve to such cool temperatures, instead traversing an 04- 
61a— >04-6Ia + — >WN7-9ha pathway, the properties of Cyg OB2 #12 are therefore difficult to understand under current evolutionary 
paradigms. Finally, we note that as with AG Car in its cool phase, despite exceeding the HD limit, the properties of Cyg OB2 #12 
imply that it lies below the Eddington limit - thus we conclude that the HD limit does not define a region of the HR diagram inherently 
inimical to the presence of massive stars. 

ON ' Key words. stars:evolution - stars:early type - stars :supergiant 
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O , 1. Introduction exploit all aspects of these stars; a problem compounded by com- 

paratively weak observational constraints on a number of key 

■ Very massive stars are of considerable astrophysical interest physical proce sses that drive massive stellar evolution, such as 
^ ; given their role in driving galactic evolution via the copious pro- the influen ce of (metallicity dependent) mass loss, the role of 

. rH . duction of ionising radiation and the deposition of chemically rotational mixing and the effect of binary interactions. In this 

^ ■ processed material and mechanical energy into the interstellar regard> a critical test of evolutionary the or y is the accurate re- 

H . medium. They are thought to be the progenitors of Type Ibc production of the properties of stars within the short lived tran- 

■ and II supernovae (SNe) and, in low metallicity environments, sitional phase betW een the main sequence (MS) and H-depleted 
gamma ray bursts (GRBs). The diverse nature of their core col- Wolf-Rayets (WRs), which is populated by a diverse 'zoo' of 
lapse SNe is also mirrored in the nature of their relativistic rem- disparate objects such as blue and yellow hypergiants (B/YHGs) 
nants which, for a given progenitor mass, may be either a neu- and luminous blue variab les (LBVs)/P-Cygni supergiants. For 
tron star or black hole depending on their pre-SNe (binary medi- example, a crucial observational finding that must be replicated 
ated?) mass loss history (e.g. Ritchie et al.|2010|>. Moreover, with is the apparent dea rth of cool evolved stars at high luminosities 
My ~ -1 1 during particular phases of their post-main sequence (the empirical Humphreys-Davidson (HD) limit; Humphreys & 
evolution, they have the potential to probe the distance to - and Davidson [1979^ in this regard the YHG/RSG populations of 
chemical composition and star formation history of - external M31 (Drout et al.|2009) and the Galactic cluster Westerlund 1 
galaxies out to distances of ~ 10 to 30 Mpc with the current and (clark et aL [20To|) already pose problems for current theoretical 
next generation of ground based telescopes (Kudritzki 2010]). models 

Unfortunately, despite considerable theoretical efforts (e.g. 
Meynet & Maeder 2005 , Heger et al. 2000) we currently lack Another manifestation of this uncertainty is an inability to 

a comprehensive theoretical framework to fully understand and incorporate the various members of the post-MS 'zoo' into a co- 
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Velocity (km/s) 

Fig.l. Comparison of three (1999-2008) and five (1994-2009) 
epochs of Ha observations for Cyg OB2 #12 (upper panel) and 
f 1 Sco (lower panel) demonstrating line profile variability at low 
projected velocities. Note the corresponding lack of variability 
in the line wings and the C ii/U 6578, 6582 doublet in the red 
wing. Due to their lower resolution, spectra of Cyg OB2 #12 
from 1992+5 are not presented, but are consistent to these spec- 
tra when all are convolved to a matching resolution. 



herent evolutionary scheme as a function of initial stellar mass. 
Various authors (e.g. Langer et al. 119941 Crowther et al. 1995 
Martins et al. 2007 2008) have proposed scenarios which have 
increasingly been informed by quantitative non-LTE model at- 
mosphere analysis of different stellar populations, latterly lo- 
cated within coeval young massive Galactic clusters. In order 
to build on this approach we present an analysis of the proper- 
ties of the BHGs of early B spectral type, which hitherto have 
escaped systematic study. 

The Ia + luminosity class was first applied to a group of four 
highly luminous (My < -8) B 'super-supergiants' within the 
LMC by Keenan d 1 97 It . which were later described as 'hyper- 
giants' by van Genderen et al. (1982). These are further dis- 
tinguished from normal BSGs by the presence of (P Cygni) 
emission in the Balmer series (cf. HD 190603; Lennon et al. 
1 19921 1. Currently, to the best of our knowledge only 16 (candi- 



date) BHGs have been identified within the Galaxy, of which 
eight are of early (<B4) and eight of late (>B5) spectral type; 
these are summarised in Table 1 . Of the early BHGs, Kaper et al. 
(2006) have already presented the results of a quantitative anal- 
ysis of BP Cm (= Wray 977; the mass donor in the High Mass 
X-ray Binary GX301-2); in this paper we present the result of 
equivalent analyses of Cyg OB2 #12, £ Sco(=HD 152236) and 
HD 190603. Foreshadowing Sect. 4, while the latter two objects 
have previously been subject to similar studies (e.g. Crowther et 
al. [2006 Searle et al. 120081 1. we employ more extensive multi- 
wavelength datasets that enable significantly more accurate de- 
terminations of physical parameters; for example the lack of 
wind contamination in the higher Balmer and Paschen lines sam- 
pled here permitting a more robust determination of the surface 
gravity. Insufficient data exist to undertake comparable modeling 
of the remaining early B hypergiants, although we are currently 
in the process of obtaining the requisite observations. 

Of these stars, Cyg OB2 #12 is of particular interest as it 
has long been recognised as one of the intrinsically brightest, 
and potentially most luminous stars in the Galaxy (=Schulte 12 
; Sharpless fTMTl Schulte [T958J, lying well above the HD limit 
and so should provide a stringent test of current evolutionary the- 
ory. Moreover, BHGs have previously been associated with the 
LBV phenomenon by various authors (e.g. Clark et al. [2005a) 
and LBVs in turn have been implicated as both critical to the 
formation of WRs and as the immediate precursors of type II 
SNe (Smith & Conti|2008ll. Finally, Cyg OB2 #12 and ^ Sco 
are generally thought to be located within the Cyg OB2 and Sco 
OB 1 associations respectively, which have both benefited from 
multiple (recent) studied and hence in principal should help in 
the assessment of both the nature of their progenitors as well as 
their placement in a post-MS evolutionary sequence. 

The paper is structured as follows. In Sect. 2 we briefly de- 
tail the new observations of the 3 BHGs analysed in this work 
(Cyg OB2 #12, £ 1 Sco and HD 190603; we utilise the dataset 
and analysis of Kaper et al. d20061 > for BP Cm), while sum- 
marising the consolidated observational datasets in Sect. 3. This 
enables us to assess the degree of short and long-term variability 
demonstrated by the program stars as well as providing the spec- 
tra and optical-radio spectral energy distributions (SEDs) for the 
non-LTE model atmosphere quantitative analysis. The results of 
this are presented in Sect. 4, discussed in Sect. 5 and conclu- 
sions presented in Sect. 6. Finally, an extensive summary of the 
(historical) observational properties of Galactic BHGs is docu- 
mented in Appendix A, the spectropolarimety of Cyg OB2 #12 
and f 1 Sco discussed in Sect. 4 is presented in Appendix B and 
Appendix C contains an analysis of the properties of the host 
cluster (NGC 6231) and association (Sco OBI) of C Sco. 

2. Observations and data reduction 

In order to accomplish the goals of the paper we have compiled 
extensive datasets for Cyg OB2 #12, C Sco and HD 190603, 
utilising both new and published observations. A full presenta- 
tion of the data for these (and other) BHGs may be found in 
Appendix A, while we briefly describe the new observations un- 
dertaken and data reduction employed below. In appendix B we 
also present the results of analysis of previously unpublished 
spectropolarimetric observations of Cyg OB2 #12 and £ l Sco; 

1 e.g. Cyg OB 2: M assey & Thomp son ( flWT) , Knodlseder (20001, 
Comeron et al. J2002b, Ha nson (20031 & Neguerue la et a l. d2008t and 
Sco OBI: R eipurth f2D08t . Sana et al. (2006a 20071 [20081 and Raboud 
et al. fl997l >. 
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Table 1. Summary of Galactic BHG properties. 



Star 


Spectral 


Co-ordinates 


Vband 


Cluster/ 




Type 


R.A. 


Decl. 


(mean) 


Association 


BP Cm 


Bl Ia + 


12 26 37.56 


-62 46 13.2 


10.68 


- 


HD 169454 


Bl Ia + 


18 25 15.19 


-13 58 42.3 


6.65 


- 


C Sco 


B1.5 Ia + 


16 53 49.73 


-42 21 43.3 


4.78 


Sco OBI 


HD 190603 


B1.5 Ia + 


20 04 36.17 


+32 13 07.0 


5.66 


- 


HD 80077 


B2.5 Ia + 


09 15 54.79 


-49 58 24.6 


7.57 


Pismis 11 


Cyg OB2#12 


B3-4 Ia + 


20 32 40.96 


+41 14 29.3 


11.47 


Cyg OB2 


Wdl-5 


WNL/B Ia + 


16 46 02.97 


-45 50 19.5 


17.49 


Wdl 


Wdl-13 


WNL/B Ia + 


16 47 06.45 


-45 50 26.0 


17.19 


Wdl 


Wdl-7 


B5 Ia + 


16 46 03.62 


-45 50 14.2 


15.57 


Wdl 


Wdl-33 


B5 Ia + 


16 47 04.12 


-45 50 48.3 


15.61 


Wdl 


HD 183143 


B7 Iae 


19 27 26.56 


+ 18 17 45.2 


6.92 




HD 199478 


B9 Iae 


20 55 49.80 


+47 25 03.6 


5.73 




HD 168625 


B8Ia + 


18 21 19.55 


-16 22 26.1 


8.44 




Wdl-42a 


B9 Ia + 


16 47 03.25 


-45 50 52.1 




Wdl 


HD 168607 


B9 Ia + 


18 21 14.89 


-16 22 31.8 






HD 160529 


B8-A9 Ia + 


17 41 59.03 


-33 30 13.7 







Note that no K-band magnitude is available for Wdl-42a, while both HD 168607 and 160529 are known large amplitude photometric variables 
(Sect. A.6). 



while these are not directly employed in the quantitative model- 
ing of these stars they serve as a useful check on the geometry 
of the cirumstellar environment of these stars. 



2.1. CygOB2#12 

New high resolution and S/N spectra of Cyg OB2 #12 were 
obtained through the William Heschel Telescope (WHT) ser- 
vice programme on 2008 July 28, employing the ISIS double- 
beam spectrograph equipped with the R1200B and R1200R grat- 
ings; a summary of the resultant wavelength ranges is given in 
Table 2. These observations were supplemented with a further 
five epochs of previous unpublished archival observations dat- 
ing from 1992-2007; these too are summarised in Table 2. All 
spectra were reduced in a consistent manner using the Starlink 
packages Figaro and Kappa, with selected regions presented in 
Figs. 1-2 and A.l. Continuum normalisation was accomplished 
via spline fitting and division, yielding uncertainties of a few 
percent at most. Likewise, the RMS on the fits to the arc spectra 
are of the order of ~0.03 pixels, resulting in a negligible error 
for wavelength callibration. 

Unfortunately the distance and reddening to Cyg OB2 #12 
preclude UV observations, so the 4000-9000A spectrum ob- 
tained in 2008 formed the primary dataset employed in this 
study. This was supplemented with the ljt/m spectrum of Conti 
& Howarth (1999), which encompasses the important Hei 
1.083jum transition, as well as the flux callibrated ISO-SWS 
spectrum presented by Whittet et al. (1 19971 1 and an archival 
Spitzer Space Telescope IRS spectrum (Houck et al. 2004). 
An optical - radio spectral energy distribution (SED) was con- 
structed from continuum flux measurements from the literature 
(Tables A.1-A.3). While these observations were not contem- 
poraneous, as demonstrated in Sect. A.l there is no evidence 
for substantial variability/evolution over the timeframe spanned 
by these observations. Hence we are confident that such an ap- 
proach is well justified and this is supported a posteriori by the 
excellent fits to the combined dataset. 



2.2. £ SCO 

A total of 18 new spectra have been utilised for the analy- 
sis of the long-term behaviour of ( l Sco. Unpublished ESO 
2.2m/FEROS spectra from 1998-9 "(Program ID 063.H-0080; 
PLKaufer) and 2005 (Program ID 075.D-0103(A); PLDufton) 
were kindly provided by Otmar Stahl (2010; priv. comm.) and 
Phillip Dufton (2010; priv. comm.) respectively. Data reduction 
was accomplished via the custom pipeline described in Stahl et 
al. ([1999 ). Continuum normalisation was accomplished by first 
dividing the spectra by the instrumental response curve, with a 
subsequent division by a spline function; as before, the errors 
in this procedure are of the order of a few percent. The spec- 
trum from 2006 February 17 was obtained by us with the ESO 
NTT/EMMI; standard reduction procedures were employed and 
are discussed in Negueruela et al. (in prep.). Finally new blue 
and red end spectra from 2009 March 06 were obtained from the 
ESO archive (Program ID 082.C-0566(A); PLBeletsky). These 
were obtained with the VLT/UVES with cross disperser gratings 
CD1, 2 and 3 and were reduced using the standard ESO pipeline. 
A full summary of these new observations is provided in Table 2 
and selected spectra are presented in Figs. 1-2 and A.l. The full 
dataset used for modeling also includes archival IUE and ISO 
spectra, as well as published photometry summarised in Tables 
Al-3. 

2.3. HD 190603 

Finally, we make use of solely archival data for the quanti- 
tative analysis of HD190603. As with the other objects, the 
sources of data used in constructing the SED are given in 
Tables Al-3. Spectroscopic data were taken from Lennon et al. 
(3950-5000A+ Ho-; [1992] ), Rivinius et al. (4050-6800A; [T997j ), 
Andrillat et al. (8390-8770A) as well as one from the IACOB 
database (3710-6800A; Sergio Simon-Diaz, priv. comm. 2011). 

3. The observational properties of BHGs 

BHGs form a spectroscopically homogeneous class of objects 
that are distinguished from BSGs by the presence of Balmer 
emission lines. However, they have been implicated in the LBV 
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Fig. 2. Comparison of synthetic spectra (yellow dashed line) 
with the optical to radio SEDs of CygOB2 #12 (top panel), 
£ 1 Sco (middle panel) and HD 190603 (bottom panel). Note the 
lack of a near-mid IR excess signalling the absence of signif- 
icant amounts (warm) circumstellar dust. Physical parameters 
employed in the modeling - such as Lboi and E(B - V) - may 
be found in Table 4 and Sects. 4.1 and 4.2. 



phenomenon (Sect. 1) and so it is of interest to compare their 
0.4-4. lyum spectra with those of bona fide LBVs in both hot and 
cool phases. Suitable comparison spectra across this wavelength 
range are provided by, amongst others, Stahl et al. (200TJ, Hillier 
et al. ST9m . Groh et al. (120071 . Clark et al. (I2003al 120091 l20TTb 
and Lenorzer et al. (2002). Significant differences are evident; 
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Fig.3. Hipparcos (Top panel:1989-1993) and NSVS (Bottom 
panel: 1999-2000) lightcurves of Cyg OB2 #12. Both panels plot 
the (non standard) instrumental magnitudes (references to these 
are provided in the text) and we caution that the scale of both 
axes differ between the two panels. 

Table 2. Summary of unpublished spectral observations. 



Date 


Telescope 


Wavelength 
Range(A) 


Resin. 




Cyg#12 






1992 July 17 


INT 


4036-4836 


4000 






5765-9586" 


4000 


1995 July 17 


WHT 


6172-6980 


4000 


1998 Aug. 09 


WHT 


3940-5737 


8000 






6366-6772 


8000 


2000 July 17 


INT 


3940-5737 


8000 






6327-6727 


8000 


2007 August 22 


WHT 


7600-8900 


10000 


2008 July 22 


WHT 


4000-4700 


4500 






4500-5270 


4500 






6450-7150 


7000 






8350-8900 


10000 




C Sco 






1998 October 7 


ESO 2.2-m 


4000-8950* 


48000 


1999 July-Aug. 


ESO 2.2-m c 


4000-8950* 


48000 


2005 April 24 


ESO 2.2m 


3800-6800 


48000 


2006 February 17 


NTT 


3933-7985 


9840 


2009 March 06 


VLT 


3060-5600 





"Note this wavelength coverage was obtained via 5 overlapping ob- 
servations. b Approximate usable wavelengths. c 14 spectra obtained be- 



tween July 17 to August 4. 



specifically the weakness of the H i and He i lines in the BHG 
spectra, as well as the absence of emission in low excitation 
metallic lines such as Fen and m that characterise the spectra 
of LBVs. 

A second observational characteristic of some LBVs is the 
presence of circumstellar dust - indeed the presence of a de- 
tached circumstellar nebula has been advanced to support the 
LBV nature of the late-B hypergiant HD 168625. Of the early 
BHGs in question, only BP Cm is known to show evidence 
for warm dust (e.g. Fig. 2; Moon et al. 2007), although with 
Miust ~ 5.2xlO~ 8 M , considerably less than is associated with 
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either LBVs or supergiant B[e] stars (Egan et al. 2002 , Clark et 
al. I2003al Kastner et al. [2006b . 

However, the defining characteristic of LBVs is their dra- 
matic spectroscopic and photometric variability over timescales 
of months-years. Motivated by reports of spectral variability 
amongst both early- and late-B HGs we undertook an exhaustive 
literature search for these objects, which we present in Appendix 
A, and briefly discuss here. As highlighted by previous authors, 
there is strong observational evidence for a physical associa- 
tion between LBVs and low-luminosity, late-B HGs such as 
HD 160529 and HD 168607; however, the link is far from proven 
for earlier, more luminous BHGs. Owing to their luminosities 
extensive spectroscopic and photometric datasets exist for such 
stars, with the majority extending back to the mid-20th century, 
for Cyg OB2 #12 and HD 190603 to the turn of last century and 
(sparse) photometric data for ( l Sco potentially many centuries 
before that (Appendix A |3 

Cyg OB2 #12 has long been suspected of being spectroscop- 
ically variable, being variously classified as B3-8 Ia + in the liter- 
ature (e.g. Table A. 4 and refs. therein). However, upon close ex- 
amination of the spectra and classification criteria employed, we 
conclude that long-term evolution of the stellar temperature be- 
tween 1954-2008 is probably absent (Appendix A. 1.2) and that 
the reports of such behaviour result from comparison of low S/N 
and resolution spectra and, crucially, the sensitivity of the com- 
monly employed He 1 4471/Mgn 4481 criterion to the properties 
of the photosphere/wind transition zone as well as photospheric 
temperature (Sect. 4.1). This conclusion is bolstered by the lack 
of secular photospheric variability dating back to the 1890s, al- 
though low level, apparently aperiodic variability does appear 
present (e.g. Gottlieb & Liller l 19781 1. 

Based on the datasets discussed in Appendix A, similar 
conclusions may be drawn for the remaining early- (BP Cm, 
HD 80077, HD 190603, HD 168454, Wdl-5 and 13) and mid- 
B HGs (Wdl-7 and 33), with the potential exception of £ l Sco. 
Spectroscopically, it is the best sampled of all the BHGs and 
shows no evidence for evolution between ~ 189 1-2009 (Table 
A. 5). Likewise, it appears to have remained photometrically sta- 
ble since at least 1949 (Table A.l), although isolated historical 
photometric observations dating back many centuries are avail- 
able in the literature, from which Sterken et al. ( 1997) infer pos- 
sible LBV like variability in the 18-19th century. 

However, while long-term variability appears rare or absent 
amongst the early-B HGs, rapid (~day to day) line profile vari- 
ability (LPV) in both wind dominated emission and photospheric 
absorption lines appears ubiquitous, persistent and well docu- 
mented (e.g. Fig. 1 and Appendix A and references therein). 
These behaviours are thought to arise from time variable wind 
structure and photospheric pulsations respectively, with the for- 
mer also potentially giving rise to variable radio continuum 
emission (e.g. Gonzalez & Canto 2008), although changes in the 
the degree of ionization of the wind are another possible cause 
of this phenomenon. 

Likewise, rapid, low amplitude photometric variability is 
present amongst all the BHGs (e.g. Fig. 3 and Appendix A). 
Such behaviour has long been recognised as being characteristic 
of luminous stars of all spectral types (the 'a Cygni' variables; 
e.g. Burki [T978l van Leeuwen et al. [19981 Clark et al. [20T0T > and 
again is typically attributed to photospheric pulsations. 



2 Unfortunately, in many cases early observations are presented in the 
literature without precise observational dates. Nevertheless, the similar- 
ity of these data to more recent observations supports the lack of long 
term variability in such cases. 



Table 3. Model Atoms included in our calculations. 



Species 


# full levels 


# super levels 


# of transitions 


Hi 


30 


30 


435 


He i 


59 


59 


592 


Hen 


5 


5 


10 


Cn 


80 


45 


648 


Cm 


38 


22 


149 


Civ* 


8 


5 


22 


Ni 


63 


25 


344 


Nil 


179 


86 


1936 


Nm 


20 


11 


61 


Oi 


75 


23 


450 


On 


102 


46 


992 


Nai 


44 


18 


345 


Mgn 


50 


37 


597 


Aln 


20 


12 


41 


Aim 


45 


17 


362 


Sin 


72 


35 


491 


Sim 


83 


47 


532 


Si iv 


28 


17 


129 


Sn 


85 


26 


745 


S in 


63 


31 


331 


Can 


32 


25 


157 


Crm 


209 


36 


2444 




Z4Z 




4 / J A 


Fen 


709 


233 


18011 


Fern 


812 


87 


21383 


Feiv 


417 


62 


8353 


Nin 


254 


33 


3354 


Nim 


336 


34 


5181 


Ni iv 


294 


56 


4174 


Total 


4160 


1137 


72827 



Not used in Cyg OB2#12. 



We conclude that the temporal behaviour of the early-B 
HGs appears entirely typical of the wider population of lumi- 
nous, early-spectral-type non-LBV stars, with the potential ex- 
ception of ( l Sco in the 18-19th centuries and a possible sudden 
Aitib ~0.4 mag 'glitch' in the lightcurve of Cyg OB2 #12 in the 
mid-1940s (Gottlieb & Liller [19781 . 

4. Quantitative modeling 

In order to more fully address the physical properties of 
Cyg OB2 #12, Sco and HD 190603 we have employed the 
non-LTE model atmosphere code CMFGEN (Hillier & Miller 
1998, 119991 1. This solves the radiative-transfer equation for a 
spherically symmetric wind in the co-moving frame under the 
constraints of radiative and statistical equilibrium. Since the 
code does not solve for the wind structure from first (physical) 
principles, a velocity structure must be chosen; we adopted a 
standard /3-type law[j. Given the relatively high wind densities 
of BHGs, the role of the hydrostatic structure and the transi- 
tion region between photosphere and wind becomes crucial to 
interpret the spectra (e.g. Sect. 4.1 and Fig. 8). Because of this, 
we used a hydrostatic density structure at depth (the 'pseudo- 
photosphere') matched to a yS-law wind, requiring the density to 
be continuous at the transition. The location of the transition is 
fixed by the adopted base velocity of the wind, v p h ot (a free pa- 
rameter in the fitting), which specifies the density through the 

3 The /3-type law is given by v(r) = v«, X (1 — b/rf, where b = 
'phot x (1 _ (v p hot/Voo) <1//?> ) and r phot and v phot are the radius and velocity 
at the connection point between the /J-law and the hydrostatic structure. 
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4107 



4123 



0.6 

4140 4320 



4347 



4360 4835 



4868 



4885 



2.5 



2.0 



1.5 



1.0 



0.5 

^ 6540 





Ha 


■ 1 ■ 




; Obs 
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1 1 

J\ \ 




Cll 






, i , 





1.0 



0.9 



0.8 



Hel 




MM 


gll 








i iv w 1 













6557 



6573 



T3 
N 

E 
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Fig. 4. Comparison of the synthetic spectrum of Cyg OB2 #12 (red dashed line) with observational data for various transitions 
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Table 4. Model parameters for the BHGsCyg OB2 #12, f 1 Sco, HD190603 & BP Cm (upper panel) and Galactic B supergiant 
comparison (lower panel; Searle et al. 2008 ) 
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Note that R, corresponds to R(t Ross = 2/3), the H/He ratio is given by number and v phot in terms of v solm d- Abundances are relative to solar values 
from Anders & Grevesse ( 19891 and have uncertainties of, typically, ~0.2dex; if we use the values from Asplund et al. ( 2006) as a reference the 
derived ratios need to be scaled by 1.49, 1.86 and 1.86 for C, N and O respectively. The results presented here assume a distance of 1.75 kpc 
for Cyg OB2 #12 (DM = 11.21). Apart from the assumed distance of 1.64 kpc for f 1 Sco (DM = 11.07, Sana et al 2006) we also provide 
corresponding stellar properties for the other distance estimate (1.99 kpc, DM = 1 1.50) assumed in previous spectroscopic studies. For HD 190603 
we display the stellar properties assuming a distance of 1.57 kpc (DM = 10.98) as well as those which would result if the object had the same M v 
as f 1 Sco. Finally a distance of 3.04kpc (Kaper et al 2006) is assumed for BP Cru 



equation of mass continuity. This approach is similar to that of, 
e.g., Santolaya-Rey et al. ( 1 19971 1, except that in our implemen- 
tation the flux-weighted mean opacity is used in preference to 
the Rosseland mean opacity Several comparisons using 'exact' 
photospheric structures from TLUSTY (Hubeny & Lanz 1995) 
showed excellent agreement with our methodQ. 

The main atomic processes and data sets are discussed in 
Hillier & Miller ( [19981 see Dessart & Hillier |20T0| for an up- 
dated description of the atomic data). The list of model atoms 
utilized in our calculations is provided in Table [3] including 
for each ion the number of super and full levels and number 
of bound-bound transitions. The CMFGEN model is then pre- 
scribed by the stellar radius, the stellar luminosity, L», the 
mass-loss rate, M, the velocity field, v(r) (defined by fl, v p hotand 
the terminal velocity, Voo), the volume filling factor, /, charac- 
terizing the clumping of the stellar wind, and elemental abun- 
dances. The interstellar reddening parameters E(B - V) and Ry 
are also obtained by fitting the model SEDs to the available pho- 
tometric data. Finally, the equatorial rotation velocity, v e sin i, 
and the atmospheric macroturbulent velocity, v macro , were esti- 
mated in Fourier space (Simon-Diaz & Herrero [20071 ) from a 
selected sample of photospheric lines. 

Given the number of free parameters required to specify a 
fit it is not possible to survey the full parameter space system- 
atically in order to establish robust error estimates; hence errors 
quoted in this paper represent the range of values for which an 
acceptable fit to the data may be obtained. Nevertheless, we dis- 
cuss error estimates for key parameters for each of the stars anal- 



4 The latest version of CMFGEN, which allows the user to com- 
pute the exact hydrostatic structure, showed full consistency with our 
method. At this stage, we prefer our approach as it allows for more 
flexibility to study the crucial transition region. 



ysed below. The validity of this technique has been demonstrated 
by calibration to stars of similar temperature and luminosity for 
which UV, optical and near-IR data were available (Najarro et 
al. fT999l Naiarro l200TT ). 

A comparison of the predicted SEDs to our data for the 
three stars is presented in Fig. [2] while selected regions of the 
synthetic spectra are overplotted on the observed 4090-1 1050A 
spectrum of CygOB2 #12 in Fig.S the UV-near IR (-1200- 
22100A) spectrum of ( x Sco in Figs. [5] & |6] and the optical 
(-4000-7000A) spectrum of HD 190603 in Fig. 7. The optical 
data shown in Fig. |6] correspond to the 2006 ESO EMMI/NTT 
run. Fits to the earlier FEROS spectra supplied by Otmar Stahl - 
which provide a more extended wavelength coverage - are shown 
in the Appendix (electronic version only) in Figs lD.ll and lD.2l 
The latter were used to obtain the final stellar properties as they 
encompass the high Balmer and Paschen lines which constrain 
the surface gravity. A summary of the results of the analyses for 
the 3 stars is presented in Table [H along with the parameters of 
the Bl Ia + star BP Cru, obtained via an identical methodology 
by Kaper et al. J2006I ). 

4.1. CygOB2#12 

For the purposes of modeling Cyg OB2 #12 we employed 
the dataset described in Sect. 3 and Appendix A. We have 
also assumed membership of CygOB2 and hence a distance, 
d~1.75 kpc (e.g. Negueruela et al. 120081 and refs. therein); we 
return to this issue below. 

Temperature, gravity and Luminosity The available optical and 
NIR spectra allowed us to make use of several ionization equilib- 
ria to estimate the effective temperature, and thus constrain the 
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Fig. 8. Importance of the wind/photosphere transition region on the emergent spectrum of Cyg OB2 #12. Left: Influence of the 
transition velocity (in units of the sound speed, v SOU nd ) on the spectral typing diagnostic lines Hei 4471 and Mgn 4481 A. Right: 
Influence on the Fe m/Fe n ionization equilibrium diagnostic lines. The Fe n 5 170 A line is severely affected by the transition velocity 
while the Fe in lines remain roughly unaltered. 



ionization structure of Cyg OB2 #12. Hence, we were able to uti- 
lize simultaneously Si in/Si n, O n/O i, N n/N i and Fe m/Fe n line 
ratios. The S m/S n equilibrium was only used for a consistency 
check, as the weak diagnostic S in 4253A line suffers from a 
poorer S/N ratio in the blue spectral region due to reddening. 
We refer the reader to Sec IA.1.2l for a detailed description of the 
available diagnostic lines. 

When analyzing the ionization equilibria, we found that 
while lines belonging to higher ionization stages react sensi- 
tively to effective temperature and to lesser extent to gravity, 
those corresponding to the lower ionization stage display a high 
sensitivity not only to temperature and gravity, but also the lo- 
cation of the transition region between photosphere and wind, 
characterised by v p hot (especially the Sin, Fen and Oi lines). 
This sensitivity is clearly illustrated in Fig. |8](right panel), where 
the choice of v p h t significantly affects the resulting strength 
of the Fen 5170A line, leaving the Fern lines basically unaf- 
fected. Thus, if a lower transition velocity is chosen, a lower ef- 
fective temperature is required to match the Fe m/Fe n ratio. As 
we will show later, this effect is responsible for our somewhat 
lower r e ff values when compared to other studies (see Sect l4.2l >. 
Interestingly, Fig. [8] (left panel) also shows how the ratio of the 
diagnostic lines Hei 4471 and Mgn 4481 A, which are used for 
spectral typing, also shows a moderate dependence on the tran- 
sition region. Thus, a model with v p h t = 0.5v sount j will yield 
equally strong He 1 447 1 and Mg n 448 1 A, while if the transition 
takes place at v p h ot = 0.05v soun d, Hei 4471 A clearly becomes 
stronger than Mg n 448 1 A. 

We note that this moderate dependence on the transition 
region becomes important in BHGs, due to the presence of a 
strong, dense wind and should be negligible in B supergiants, 
where the absorption lines will form in deeper photospheric lay- 
ers that are largely unaffected by the photosphere- wind transition 



region. Therefore, despite the large number of diagnostics, our 
final estimate of the effective temperature, (!T e ff = 13.7kK, see 
Table [4} is subject to a moderate uncertainty. We find +800K 
and -500K as upper and lower error bounds. 

Once the effective temperature was obtained, and assuming 
a distance of d~1.75 kpc, we proceeded to fit the observed SED 
of Cyg OB2 #12 from the optical through radio (see Figf2]i and 
hence derived the reddening, stellar radius and, therefore, the 
stellar luminosity. We found E(B - V) = 3.84 and a reddening 
parameter Ry = 2.65, corresponding to Ay = 10.18. This value 
shows excellent agreement with the Ay - 10.20 found by Torres 
et al. 119911 although they obtained a higher total to selective 
extinction parameter (Ry = 3.04) 

From these values we obtained a stellar radius of 246/? Q and 
a final luminosity of 1.9 x 10 6 L Q . The temperature and resul- 
tant luminosity are broadly comparable to previous qualitative 
studies (e.g. Massey & Thompson 1991] Hanson[2003), confirm- 
ing that it is an extraordinarily luminous, but comparatively cool 
BHG. 

The Paschen lines in the I Band provide the best constraints 
for the surface gravity, especially the run of the line overlap 
among the higher members (Fig. 0J. Compared to the classi- 
cal optical diagnostic lines - Hy and HS - the higher Paschen 
series lines are significantly less affected by the stellar wind 
(and consequently clumping) and the wind/photosphere transi- 
tion region. Nevertheless, the influence of the latter once again 
translates into a larger uncertainty in the (lower) error. Thus, we 
find \ogg = 1 -70*qj5, corresponding to a spectroscopic mass of 
M=\Wf x M Q . 

Wind properties and clumping Unfortunately, the moderate 
reddening affecting Cyg OB2 #12 prevents us from securing UV 
observations from which one might derive firm v m estimates. 
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Fig. 5. Comparison of the synthetic (red dashed line) with the 
observational (black line) UV spectrum of ( l Sco. Note that the 
narrow lines that are not replicated in the model are interstellar 
in origin. 



Furthermore, as with Sco and HD 190603, Cyg OB2 #12's 
stellar wind is not sufficiently dense that alternative v M diagnos- 



tic lines such as Ha or Hei 10830A can reach their full poten- 
tiallity and unambiguously yield the wind terminal velocity. 

Highlighting this difficulty, Souza & Lutz ( 1980) suggested 
a terminal velocity of 1400km s , based on the presence of a 
blueshifted absorption feature in the Ha profile; a value which 
was subsequently used to obtain the mass-loss rate from radio 
(Abbott et al. |1981) and IR measurements (Leitherer et al. |1982) . 
Subsequently, Klochkova & Chentsov (2004) revised the esti- 
mate of the wind's terminal velocity significantly downwards. 
Based on higher resolution Ha observations they identified 
strong electron scattering wings extending up to 1000km s _1 , as 
well as blueshifted absorption up to ~ 160km s _1 which they at- 
tributed to the wind's true v x . However our models indicate that 
this blueshifted feature does not directly reflect but rather re- 
sults from the run of the density and ionization structures within 
the wind which shapes Ha. 

This is shown in Fig. 9, where models with v m ranging 
from 400 to 1000km s _1 yield an absorption feature around 100- 
200km s _I to the blue of Ha in the resultant synthetic spectra. 
Moreover, values of significantly below 400km s appear 
too low to reproduce the high velocity line emission while, de- 
pending on the normalization errors of the Ha profile, values 
between 400 and 1000km s _1 are potentially consistent. We also 
found that the blue absorption component of the Hei 10830A 
line is likewise not able to distinguish between Vco = 150kms _1 
and Voo = lOOOkmi -1 . We note that if the wind density of 
Cyg OB2 #12 was a factor of two higher, the blue absorption 
component would develop up to the corresponding value of v M . 

Nevertheless, an upper limit on may be obtained by 
means of the Bra/He i 4.048/im complex. This is apparent in 
Fig-EJ where we see that ISO observations of this feature clearly 
resolve both components. Our models indicate that if v ro is above 
~500km s _1 , both components are blended, while strong emis- 
sion bluewards of 400km s~'is present, which is not seen in the 
data. On their own, these observations lack sufficient S/N and 
spectral resolution to accurately constrain v x , but in conjunction 
with Ha they jointly provide stronger constraints. Following this 
approach, we adopt v M = 400kms~ 1 for the remainder of the pa- 
per, while recognising that values between 300 and 1000km s _1 
cannot formally be discarded at present. Regarding this, we sug- 
gest that ground based high resolution spectra of Bra and with 
sufficient S/N to trace the line wings may constitute the best ob- 
servational constraint available to estimate v M . 
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Likewise, we make use of the shape of Ha and the 
Bra/He i 4.048/mi complex to estimate /3, the parameter charac- 
terizing the velocity law. We obtain [} — 3 indicating a relatively 
flat velocity field. Values of /3 below 2 or above 4.5 can be ruled 
out from the line fits. 

We are, however, unable to fully reproduce the blue shoulder 
of the Ha line, noting that this is the case irrespective of the ter- 
minal wind velocity adopted in the modeling. In this regard we 
highlight comparable discrepancies between the synthetic and 
observed Ha profiles for Sco and HD 190603 (cf. an inability 
to reproduce the P Cygni absorption features in these stars; Figs. 
6 and 7) while in many cases the model fits to BSGs presented 
by both Crowther et al. (120061) and Searle et al. (120081) also suffer 
similarly. 

Regarding Cyg OB2 #12, the lack of changes in polarisation 
through the line (Appendix B) strongly argues against any large 
scale wind asymmetries that might have been supposed to ex- 
plain this feature. Hence we feel confident in the application of 
CMFGEN - which adopts spherical symmetry - to this and other 
stars in this study, while noting that this discrepancy might indi- 
cate a shortfall in the physics employed. Subject to this - and the 
preceding - caveats, we emphasise that the presence of this dis- 
agreement between model and observations does not affect our 
ability to determine the terminal velocity of the wind following 
the methodology described abov^|. 

Finally, we were able to derive a mass-loss rate of M- 
3.0 x 10~ 6 M o yr~' (Table 4). Of course this value is bound to 
the adopted v M and clumping law. Thus, if Voo values of 300, 700 
and 1000km s~'were adopted, fits of similar quality would be 
obtained for M~ 2.5, 4.0 and 6.5 x 10~ 6 Af o yr _1 respectively. 

The main observational constraints which set the run of the 
clumping law are the Ha - and Bra emission components and 
the IR and submillimeter + radio continuum. An onset of the 
clumping at relatively high velocity ( ~ 200km s -1 ) is required 
to avoid strong emission in the above lines and too great an ex- 
cess in the IR continuum. On the other hand, we find that a final 
clumping value of 0.04 is required to reproduce the submillime- 
ter and radio continuum (see Fig|2]i. 

Our derived mass loss rate presents one of the major differ- 
ences with respect to previous works. The significantly lower 
clumped mass loss rate of M~ 3.0 x 10 _6 M o yr _1 corresponds 
to an unc lumped (M// 5 ) rate of ~ 1.5 x 10~ 5 M G yr _1 ; a factor 
of ~3 lower than that reported by Leitherer et al. (1982), due 
to the considerably smaller wind terminal velocity adopted 
(400 versus 1400kms~'). The effects of the high degree of wind 
clumping found for Cyg OB2 #12 mirror recent findings for 
other massive (evolved) stars (e.g. Najarro et al. 2009, Groh et 
al. l2009all2009bl) . 

Nevertheless, such a clumping corrected value is in excess 
of that found for the less luminous and massive B supergiants 
studied by Crowther et al. (120061 ) and Searle et al. (2008j). For 
wind velocities v M < 400 kms _I the clumping corrected mass 



5 For completeness we also explored the possibility that the discrep- 
ancy between observations and model could be due to the presence 
of (spatially unresolved) blueshifted nebular emission. Fitting an ad- 
ditional simple gaussian profile to the blue emission 'shoulder' in the 
Hff and ji lines, scaled assuming typical nebular line intensity ratios 
(Ha/yfl ~2.9 and H/?/y ~2.2), resulted in significantly improved line fits. 
However, we caution that the adoption of such a solution appears pre- 
mature at present as no other expected nebular emission lines such as 
[N a], [S n] or [O i] are present in the spectrum of Cyg OB2 #12; further 
adaptive optics or coronographic observations to search for compact 
nebular emission would be of value to determine if such an approach is 
physically well motivated. 
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Fig. 10. Montage of selected He i lines in the spectrum of Cyg 
OB2 #12 with synthetic spectra computed utilising different 
He/H ratios overplotted. 



loss rate approaches those of known LBVs (Fig. 11; Sect. 5), 
although the terminal velocity is comparable to BSGs of similar 
temperature. Adopting v^, ~ 1000 kms -1 leads to a mass loss 
rate comparable to those of the LBVs, but such a wind velocity 
is significantly greater than those determined for both LBVs and 
BSGs of equivalent spectral type. With respect to this, the lack of 
the low excitation metallic emission lines that characterise cool 
phase LBV spectra (Sect. 3) is a result of the extremely large 
radius found for Cyg OB2 #12. 

Abundances Fortunately, unlike stars such as HDE 316285 
(Hillier et al. [T998i l. Cyg OB2 #12 does not suffer from a r eff /He- 
abundance degeneracy. Since, as previously described, T e g can 
be accurately determined from e.g. the Si m/Si n or O n/O i equi- 
libria, we may therefore reasonably constrain the He/H ratio. In 
Fig. 10 we overplot synthetic spectra constructed with a wide 
range of He/H abundances on selected He i transitions; note that 
as expected the best constraints are provided by the weaker, non 
saturated, transtions. This reveals one of the fundamental results 
of our study - that Cyg OB2 #12 demonstrates a solar H/He=10 
ratio (by number). Unlike the rest of BHGs analysed in this paper 
(see Table |H, which clearly demonstrate surface helium enrich- 
ment, the He i lines of Cyg OB2 #12 are best reproduced assum- 
ing no helium enrichment at all (noting that previous studies of 
BSGs adopted H/He=5 for all stars considered; Crowther et al. 
2006, Searle et al. 2008). Furthermore, even better fits to some 
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Fig. 11. Comparison of the wind properties of the BHGs to field 
BSGs (plotted according to spectral type), LBVs of tempera- 
ture selected to match the BSGs (>15kK) and those of com- 
parable luminosity to Cyg OB2 #12 (log (L/L Q ) > 6.0). Data 
from Crowth er et al. ( 120061, Se arle et al. (l20"0"8l l , Najarro et al . 
d 1999120091 ). Clark et al. (120091 and Groh et al. d2009al20 09b). 
Inclusion of later BSGs continues the trend to lower mass loss 
rates and wind velocities (e.g. Markova et al. 120081 1. Due to the 
uncertainty in the wind velocity for Cyg OB2 #12, we present 3 
pairs of possible values for Voo = 300, 400 and 1000 kms -1 (Sect. 
4.1). 
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Fig. 12. HR diagram with the positions of the BHGs (red dia- 
monds) plotted with respect to the evolutionary tracks of Meynet 
& Maeder (120051 ) and the empirical HD limit. The position 
of HD 80077 assumes membership of Pismis 11 and hence 
a distance of 3.5 kpc (Marco & Negueruela 2009). The po- 
sitions of other evolved stellar populations are also indicated: 
black squares - WN7-9h, 04-6 la and 04-6 If + stars within the 
Arches; purple circles - galactic field BSGs; blue triangles - (can- 
didate)LBVs and WNL stars; green stars - YHGs. Data from 
Martin s et al. d2007l 120081 Crowther et al. (120061 Searle et al. 
(2008) and Clark et al. (in prep.). 



of the strong He i lines in the R- (He i 6678 A), /- (He i 8581 and 
8845 A) and Z-bands (Hei 10830A) are obtained if we assume 
He to be underabundant. Moreover, we find that H/He~8.0 ap- 
pears to provide a robust lower limit to the abundance ratio, with 
lower values resulting in unacceptably poor fits to the data (Fig. 
10). 

We note that given the current stellar temperature and lumi- 
nosity, the lack of He enrichment at the stellar surfaces is at odds 
with the predictions from evolutionary models and challenges 
present theory of stellar evolution of very massive stars. We dis- 
cuss this somewhat unexpected finding, and the implications for 
the evolutionary state of Cyg OB2 #12 in Sect. 5. From Table|4] 
we see indications of CNO processing (N enhancement and CO 
depletion) in Cyg OB2 #12, though to a lesser extent than for 
( x Sco and HD 190603. This result is consistent with the derived 
H/He ratios. 

4.1 .1 . The unexpected luminosity of Cyg OB2 #12 

A critical finding of this analysis is the rather extreme luminos- 
ity of Cyg OB2 #12 (e.g. Fig. 12 and Sect. 5); are we overesti- 
mating this? An obvious explanation is that a chance alignment 
with its 'host' association leads us to overestimate its distance. 
If this were the case, the large degree of reddening demonstrated 
by Cyg OB2 #12 would either be due to a (narrow) sightline of 
anomalously high interstellar extinction towards it and only it - 



but which nonetheless yielded Diffuse Interstellar Bands (DIBs) 
comparable to those found for other bona fide association mem- 
bers (e.g. Hanson 2003 ) - or would be circumstellar in origin 
(e.g. Massey & Thompson 1991), thus permitting a more nor- 
mal run of interstellar reddening for a lower distance. However, 
our model fits to the expanded SED show no evidence for an 
IR excess at wavelengths shorter than 30/im; suggesting that any 
circumstellar component would have to be significantly cooler 
than observed around known LBVs, where pronouned emission 
from warm dust is detected at <25//m (e.g. Egan et al. 2002 
Clark et al.[2003a]). 

A second alternative is that it is a multiple system. Given 
the recent suggestion of a high binary fraction within CygOB2 
(Kiminki et al. 120091 ) and the hard, luminous X-ray emission 
from Cyg OB2 #12 (Albecete Columbo et al. 120071 Rauw et al. 
2005 ) this would appear to be physically well motivated. To re- 
duce the luminosity of Cyg OB2 #12 such that it was marginally 
consistent with current theoretical predictions (Sect. 5) would 
require a companion(s) of equal luminosity. However clear spec- 
troscopic signatures of binarity in terms of (i) radial velocity 
shifts, (ii) the presence of double lines and/or (iii) the dilution of 
spectral features appear absent- Moreover, the excellent model 
fit to the observed spectral energy distribution also excludes a 

6 In particular we note that no (anti-phased) radial velocity shifts in- 
dicative of binarity are present in the additional emission components 
of the Ha, /3 and y lines. 



13 



J. S. Clark et al.: On the nature of the galactic early-B hypergiants 




Fig. 13. HR diagram for NGC6231 (black symbols) and the 
wider Sco OBI association (blue symbols) with the non rotating 
isochrones of Schaller et al. (1992, solid lines) and the rotating 
isochrones of (Meynet & Maeder 2005 ; dashed line) overplotted. 
Open and filled symbols represent main sequence and evolved 
stars, while diamonds indicate binarity. For HD 152248 we plot 
the position of both the individual components and also that of 
the integrated system. See Appendix B. for details of contrac- 
tion. 



companion of comparable luminosity but very different temper- 
ature. 

While not fatal to a binary scenario, these observational con- 
straints would require any putative binary companion of com- 
parable luminosity to have an identical spectral type - to avoid 
its spectral signature being visible - and for the binary to either 
be very wide or seen face on - to avoid RV shifts in the spec- 
trum. However, the temperature of the hard X-ray component 
(kT = 1 .8keV) implies, via the strong shock condition, that the 
companion would have to have a wind velocity of ~1300kms~' - 
inconsistent with the emission arising in the wind collision zone 
produced by such a 'twin' and necessitating the presence of a 
third, unseen star to yield the high energy emission. While we 
cannot exclude such a 'finely tuned' hierarchical system, we find 
no compelling evidence for one either. 

Lastly, one might suppose that Cyg OB2 #12 were cur- 
rently undergoing a long-term'LBV eruption' resulting in an in- 
creased luminosity. However, while the physical properties of 
such events are ill constrained we find the current physical prop- 
erties to be distinct from those of quiescent LBVs (e.g. Sect. 
5), while mass loss rates of >> 10 _5 M Q yr _1 are typically in- 
ferred for such stars in outburst (Clark et al. 120091 1. Therefore, 
while we may not formally exclude any of these possibilities, 
we find no compelling observational evidence to support them 
either and hence consider it likely that the primary BHG within 
Cyg OB2 #12 is indeed unexpectedly luminous in comparison to 
theoretical predictions. Moreover, even if its luminosity were re- 
duced via one of these scenarios such that it was consistent with 



Fig. 14. HR diagram for Cyg OB2 with the non rotating 
isochrones of Schaller et al. ( 1992, solid lines) and the rotating 
isochrone of (Meynet & Maeder l2005l dashed line) overplotted. 
Data from Negueruela et al. (2008) and Hanson et al. (2003); the 
reader is refered to these studies for the positions of individual 
stars, which have been omitted here for reasons of clarity. 



current evolutionary tracks, we note that the revised mass loss 
rate (Moc d 15 ) and (unchanged) chemistry would still remain 
discrepant with respect to field BSGs and BHGs. 

4.2. C Sco and HD 190603 

For the analysis of [} Sco we adopted a distance d ~ 
1.64 kpc, based upon the revised estimate for NGC6231/Sco 
OBI (Appendix C); for convenience we also show results for 
the slightly larger distance of d ~ 1.9 kpc adopted in previ- 
ous works to enable a direct comparison(Crowther et al. 2006). 
Unfortunately, HD 190603 is not associated with a host clus- 
ter making the distance significantly more uncertain. For con- 
sistency to prior studies we adopt an identical distance (d ~ 
1.57 kpc; Crowther et al. 120061 1. while also presenting results 
under the assumption that it has an identical intrinsic visual mag- 
nitude to £ l Sco. 

As in the case of Cyg OB2 #12, several ionization equilibria 
could be used to determine the stellar temperature. Interestingly 
we could use three ionization stages of silicon (Si iv/Si in/Si n) 
for both objects. In the case of HD190603, our hottest object, 
the Fe n and N i lines are not detected so we used them in our 
models as indicators of a lower limit for T e ff. On the other hand, 
He ii 4686 A was weakly detected on HD 190603 allowing to use 
as well the He n/He i ionization criterium. 

Since UV spectra are available for both objects, the termi- 
nal velocity, Voo, could be accurately determined for both objects 
(e.g. Fig. 5). Given that we are fitting non-simultaneous obser- 
vations from the UV to the IR we might expect to find some 
transitions in this wavelength range which are not well fitted 
by our synthetic spectra. This appears to be the case for C rv 
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1548- 1551 A, which is significantly underestimated by our mod- 
els. Moreover these, as well as the Si iv 1394-1403A UV lines, 
are very sensitive to X-rays in this parameter domain; thus they 
should be regarded with caution. 

The main stellar properties of both objects are presented in 
Table [4] Uncertainties in stellar parameters are similar to those 
derived for Cyg OB2 #12, except for the effective temperature, 
which is slightly better constrained AT e g = +500K and the ter- 
minal velocity with Avoo = ±50kms _1 . Errors on the H/He ratio 
are of the order of ~20%. Interstellar reddening was also deter- 
mined; for f Sco we found E(B — V) = 0.66 and Ry = 3.3 and 
forHD 190603 E(B - V) = 0.70 and fly = 3.1. 

We found both stars to be less extreme than Cyg OB2 #12 
in terms of luminosity, although £ ] Sco is still significantly in 
excess of the range of luminosities spanned by Galactic 09-B5 
Ia/b supergiants (Searle et al. 2008, Crowther et al. 2006), even 
with the downwards revision of the distance. HD 190603 lies 
in the upper reaches of this distribution, but we caution that this 
could be subject to revision given the uncertainties in its distance 
estimate. As expected from their earlier spectral types, they are 
hotter than Cyg OB2 #12 , although the temperatures we find are 
slightly lower (<lkK) than previous estimates for these objects. 
These lower temperature estimates are due in part to the effect 
of the wind transition region (see Sect. 14.11) and to the use of 
all available ionization equilibria criteria. Following from these 
results, spectroscopic mass estimates are significantly lower than 
that of Cyg OB2 #12 for both stars. 

Terminal wind velocities are similar to Cyg OB2 #12, but 
lie at the lower range of values found for BSGs of comparable 
spectral type; in contrast the clumping corrected mass loss rate 
for HD 190603 lies at the upper range of those found for Bl-2 
Ia/b supergiants while that of ( l Sco is considerably in excess 
of this range (Fig. 1 1). These compare favourably with previous 
estimates for these stars (Searle et al. 2008, Crowther et al. 2006 
Markova & Puls l2008l ); although as with Cyg OB2 #12, the high 
degree of clumping required to fit the spectra provide a signif- 
icant downwards revision of absolute mass loss rates between 
this and prior studies. Values of logg are also lower than previ- 
ous estimates of both stars, due to the greater accuracy afforded 
us by the inclusion of the higher Balmer and Paschen lines in the 
spectroscopic datasets. 

Finally, in contrast to Cyg OB2 #12, but in line with esti- 
mates for field Galactic BSGs, we find a moderate H-depletion 
(or He-enrichment; H/He~5 by number) for both stars as well as 
evidence for stronger CNO processing (higher N). 

5. Discussion 

While the galactic BHG population spans over an order of mag- 
nitude in luminosity (Fig. 12), when combined with previous 
modeling results for BP Cm (Table 4; Kaper et al. 2006}, our 
analyses show the early B 1-4 Ia + hypergiants to be, as expected, 
rather luminous evolved stars, with HD 80077 and CygOB2 #12 
appearing to lie above the empirical HD limit (Fig. 12; Marco 
& Negueruela 2009). In comparison to Galactic BSGs of the 
same spectral subtypes they are overluminous and support sig- 
nificantly higher mass loss rates, with wind terminal velocities 
at the lower end of the range populated by the BSGs. 

5. 1 . The evolutionary status of BHGs 

Subject to these observational constraints how may single (early) 
BHGs be understood in evolutionary terms? Several authors 



have presented evolutionary sequences for massive (> 3OM ) 
stars, including Langer et al. (I1994I ). Crowther et al. dl995l ). 
Meynet & Maeder (2005} and Martins et al. (I2007ll2008l . After 
the initial O-type main-sequence phase, the recent empirical 
evolutionary scheme of the latter authors may be summarized as: 

• ~ 60 - 12OM : O la -> OIaf + -> WNL + abs -> WN7 

• ~ 30 - 6OM : O/B la -> WN9-llh =± LBV -> WN8 -> 
WN/C 

While not explicit in the above scheme, the properties of 
LBVs such as AG C ar (Groh et al. I2011I ) and the Pistol star 
(Najarro et al. I2009I ) strongly imply that even very massive 
Mi n i > 6OM stars may experience an LBV phase. Moreover, it 
is likely that a further subdivision will occur for M; n i < 4OM 
to accommodate the high luminosity (log(L/L Q ) ~ 5.8) cool 
hypergiants present in clusters such as Westerlund 1 (Wd 1; 
Clark etal.[2005b]): 

• ~ 30 - 40M G : O ^ OB la ^ cool Ia + ^ WN ^ WC 

noting that the evolution of spectral classification through the 
red loop and the final state reached prior to SN is particularly 
uncertain. 

Given that the physical properties of the BHGs appear more 
extreme than the field BSGs and that they appear co-located with 
the LBVs on the HR diagram (Fig. 12) a close physical associ- 
ation is suggested and we propose that the early BHGs - with 
the notable exception of CygOB2 #12 - are the immediate evo- 
lutionary progenitors of LBVs for stars of ~ 40 - 6OM . Thus, 
they would form lower mass analogues of the 04-6 Iaf + stars in 
the Arches, which Najarro et al. (120041 ) and Martins et al. |2008 ) 
demonstrate are intermediate between the 04-6 la and WN7-9ha 
stars both in terms of the degree of (increasing) chemical evolu- 
tion and also wind properties, with advancing evolutionary state 
leading to a reduction in terminal wind velocities coupled with 
an increase in mass loss rate. 

The BHGs appear to follow an analogous evolutionary tem- 
plate, with wind properties approaching those of the LBVs (Fig. 
11). Moreover, the BHGs also appear less chemically evolved 
than (candidate) LBVfl again suggesting that they occupy an 
earlier evolutionary phase. Such an hypothesis appears bolstered 
by the apparent lack of long-term(secular) variability demon- 
strated by the early-B HGs (Sect. 3), although we are aware that 
bona fide LBVs such as P Cygni may also encounter long peri- 
ods of quiescence. 

We emphasise that this evolutionary scheme is appropriate 
for single stars, but it also appears possible that early-B HGs 
may form via binary interaction and, through such a channel, 
potentially from lower mass stars. Wellstein & Langer (1999) 
present the results of simulations that show the properties of 
the known binary BP Cru are consistent with quasi-conservative 
mass transfer in a 26M Q + 25M progenitor system. Likewise, 
the BHG/WNVL transitional object Wdl-13 is the primary of a 
23.2M + 35.4M 9.27 day period binary (Ritchie et al. l20T0l ). 
This configuration is too compact for the primary to have passed 
through a BSG phase, and hence presumably formed via an 
episode of binary mediated mass loss (Sect. A.5), with a likely 



7 HD 316285 (H/He~1.5), AG Car (H/He~2.3), W243 (H/He~5), 
P Cygni (H/He~3.3), FMM3 62 (H/He ~2.8) a nd the Pisto l Star 
(H/He~1.5; Hillier et al. [19981 Groh et al. |2009a| Ritchie et al. |2009l 
Najarro et al. l2001ll2009l respectively). 
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progenitor mass of ~ 35M . In this respect we highlight that 
Wdl-13 appears spectroscopically distinct from the BHGs con- 
sidered here, showing pronounced (variable) emission lines in 
the Paschen series (Fig. A. 2 and Ritchie et al. 120101 . 

Finally, is it possible to place the later (lower luminosity) 
BHGs into comparable sequences? The consistency of the pop- 
ulation of Bl-3 la and cooler B5-9 Ia + stars within Wd 1 with 
theoretical evolutionary isochrones (Negueruela et al. 2010) sug- 
gests an extension of the paradigm that BHGs are the more phys- 
ically extreme direct descendants of hotter O/BSGs to lower 
progenitor masses (~ 35M ). However, the placement of these 
BHGs in a pre- or post-RSG/LBV phase is currently uncertain 
in the absence of tailored abundance analyses - ESO/VLT ob- 
servations currently underway will allow this to be directly ad- 
dressed. Finally, following previous suggestions in the literature 
(Sect. A. 5) - and again in the absence of quantitative analyses - 
it appears likely that lower luminosity (log(L/L ) < 5.5) BHGs 
are post-RSG objects possibly encountering an LBV phase. 

5.2. Constraints from the host stellar aggregate 

Do other observations support such an evolutionary scheme? 
Following the above discussion, the stellar populations of both 
the Arches and Wd 1 are clearly consistent. However we would 
also predict the presence of early BHGs within the Quintuplet 
cluster, which is likely to be intermediate in age between these 
clusters (Figer et al. I1999I I. While neither these authors nor 
Liermann et al. (2009) report the presence of BHGs, we note 
that e.g. LH096, 100 and 146 - which are classified as 06- 
8 If by Liermann et al. - have K band spectra identical to 
that of £ l Sco (Fig. 6). Simply applying the bolometric cor- 
rection derived from our analysis of ( l Sco yields an estimated 
log(L/L ) ~ 6.0 - 6.3 for these stars; directly comparable to the 
Pistol Star and FMM362 (Najarro et al. |2009l > and hence consis- 
tent with a pre-LBV BHG identification. Further high S/N and 
resolution observations of these stars would permit quantitative 
testing of this assertion by virtue of an accurate determination of 
temperature, luminosity and wind properties. 

Encouragingly, stars with a similar K band spectral morphol- 
ogy are also present in the Galactic Centre cluster, where they 
too are found to be of higher temperature and possess faster, 
lower mass loss rate winds than the more evolved WNL - and 
candidate LBV - stars present (Najarro et al. 1 19971 Martins et al. 
120071 

Importantly, are the properties of £ l Sco and Cyg OB2 #12 
and their host stellar aggregates consistent with this evolution- 
ary hypothesis? Unfortunately, observations of Cyg OB2 and 
NGC6231/Sco OBI raise the possibility of the stars in either 
association being non-coeval (Negueruela et al. 2008; Appendix 
C). Nevertheless, the physical properties of £ l Sco appear con- 
sistent with the population of (initially) rapidly rotating, 5 Myr- 
old stars present within NGC623 1 (Fig. 13 and Appendix B; not- 
ing that the low v e sin i of ^ Sco may be understood as a com- 
bination of spindown and/or an unfavourable inclination) or a 
somewhat younger (~ 3.2 -4 Myr) population of less rapidly ro- 
tating objects that may also be present (Appendix B). This would 
imply a progenitor mass of ~ 40 - 6OM for Sco, depending 
on age and initial rotation, broadly consistent with the current 
spectroscopic mass of ~ 36M . 

Finally, we address Cyg OB2 #12. As highlighted previously, 
this star appears difficult to accomodate in the current evolution- 
ary scheme due to the combination of high luminosity and low 
temperature displacing it from any theoretical isochrones appli- 
cable to Cyg OB2 (Fig. 14; Negueruela et al. 120081 . While its 



co-location in the HR diagram with bona-fi.de LBVs such as the 
Pistol Star, FMM 362 and AFGL 2298 (Fig. 12) suggests a sim- 
ilar nature, the long-term stability combined with its relatively 
unevolved chemsitry - and critically solar H/He ratio - appear to 
exclude such an identification. We are thus currently unable to 
place Cyg OB2 #12 into a coherent evolutionary scheme, noting 
that with a current spectroscopic mass of ~ 11OM - signifi- 
cantly in excess of both the Pistol Star and FMM 362 (Najarro et 
al. 120091 1 - one might expect it to evolve through a much hotter 
04-6I/Iaf + /WNLha evolutionary sequence. 

In this respect, the results of a quantative analysis of 
HD 80077 would be of particular interest, given its apparent sim- 
ilarly extreme luminosity (Marco & Negueruela 2009). Such an 
analysis has been performed for the M33 Bl Ia + star [HS80] 
1 10A by Urbaneja et al. ( 2005) which yielded stellar parameters 
that, like Cyg OB2 #12, both place it above the HD-limit and ap- 
pear to distinguish it from known LBVfl This finding appears to 
suggest that there might indeed be an hitherto unidentified (but 
rarely traversed?) evolutionary pathway that leads to the produc- 
tion of very luminous but comparately cool BHGs at very high 
stellar masses. 



6. Conclusions 

In order to determine the physical properties and hence evolu- 
tionary state of Galactic early-B hypergiants, we have under- 
taken new quantitative analyses of CygOB2 #12, ^ Sco and 
HD 190603 and employ previous model results for BP Cm. 
Synthetic spectra and SEDs were calculated and compared to 
comprehensive UV/optical-radio spectroscopic and photometric 
datasets compiled from the literature and supplemented by new 
and previously unpublished observations. Building on this effort, 
we also constructed exhaustive spectroscopic and photometric 
histories for all Galactic BHGs in order to test the assertion that 
they are physical identifiable with LBVs. 

Turning first to stellar variability and both the early- and late- 
B HGs demonstrated the rapid (~days) spectroscopic and pho- 
tometric variability that is symptomatic of the aspherical wind 
substructure and photospheric pulsations that characterises lumi- 
nous early stars. While low-luminosity, late-B HGs have long be 
known to undergo LBV excursions (e.g. HD 160529) we found 
no evidence for such secular variability amongst the early, high 
luminosity BHGs, with the possible exception of £ l Sco in the 
sparse photometric data obtained prior to the 2Q' h Century. In 
this regard it is interesting that both Cyg OB2 #12 and HD 80077 
appear to lie above the empirical HD limit. 

Model results reveal that the early BHGs have physical prop- 
erties - luminosity, wind terminal velocity, mass loss rate and 
chemical abundances - intermediate between Galactic BSGs and 
LBVs of comparable temperature, suggesting they provide a link 
between both evolutionary phases for stars of inital mass be- 
tween ~ 40 - 6OM . In this respect, they would play a simi- 
lar role to the more massive and luminous 04-6Iaf + objects for 
stars in the ~ 60 - 12OM range (Najarro et al. 120041 Martins 
et al. I20081 >. The simultaneous presence of populations of early 
Bl-3 la supergiants and mid-late B5-9 Ia + BHGs within Wd 1 
suggest that BHGs are potentially the immediate descendents of 
stars with initial masses as low as 3O-35M . However, the pres- 
ence of early BHGs within this cluster indicates that a parallel 
channel may also lead to this phase via close binary evolution. 



8 log(L/L G ) ~ 6.5, r eft ~ 21kK, v„ ~ 750kms-» and M= 
l.ttexlO^Moyr 1 . 
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Nevertheless, Wd 1 and the Quintuplet provide direct ob- 
servational tests of the evolutionary sequence proposed; in the 
former cluster the BHGs should be more chemically evolved 
than the BSG population but less than the RSGs that are present, 
while in the latter cluster one would expect to find a population 
of early BHGs; observations aimed at verifying both hypotheses 
will be undertaken later this year. 

The apparent physical association of Sco with 
NGC 6231/Sco OBI also allows us to test this scenario. 
Encouragingly, we find that the physical properties derived from 
our analysis are consistent with the observed stellar population 
of this region and in turn with theoretical isochrones for stars of 
between 3.2-5 Myr age - further discrimination is difficult in the 
absence of an absolute rather than projected rotational velocity. 

Conversely, the combination of extreme luminosity and cool 
temperature of CygOB2 #12 is inconsistent with theoretical 
isochrones, assuming membership of Cyg OB2. Given its co- 
location with a handful of bona fide LBVs above the HD limit it 
might be supposed to be a similar object, but a combination of 
(i) a lack of the characteristic LBV variability, (ii) extreme cur- 
rent stellar mass and (iii) lack of chemical evolution differentiate 
it from these. We are therefore currently unable to place this star 
in a consistent evolutionary scheme, noting that Negueruela et al. 
(2008 ) also find that a number of the hottest and most luminous 
stars within CygOB2 are difficult to accommodate in a forma- 
tion scenario for the association where the most recent epoch of 
star formation occurred ~2.5 Myr ago. 

Nevertheless, we note that the stability of CygOB2 #12 for 
the last century would appear to indicate that the HD limit does 
not appear to define a region of the HR diagram utterly inimi- 
cal to the presence of massive stars, but a combination of (pul- 
sational?) instabilities and extreme mass loss rates presumably 
prevent stars residing above it for large fractions of their life- 
time. Indeed, from Fig. 12 we may see that in addition to highly 
luminous BHGs such as Cyg OB2 #12 the empirical HD limit 
appears to be delineated by LBVs and other closely related (vari- 
able) stars such as the cool hypergiants and the WNVLs, further 
suggesting that it should not be regarded as a firm barrier to stel- 
lar stability, nor the location at which the onset of instabilities 
occurs. In this regard the Bl Ia + star [HS80] 1 10A and the cool 
F-hypergiant B324 (both located in M33) are of considerable in- 
terest; like Cyg OB2#12, the combination of temperatures and 
luminosities they demonstrate are not replicated by current evo- 
lutionary tracks and place both stars well above the HD limit 
(Urbaneja et al. 2005, Clark et al. 2012). Quantitative modeling 
of the Galactic BHG HD 80077 would therefore be of consider- 
able interest to see whether this too violated the HD limit. 

Moreover, despite residing above the HD limit, we find an 
Eddington parameter, TEdd, of only ~ 0.41 (and Fi-otai ~ 0.56 at 
R(t Ross = 2/3) for Cyg OB2 #12; below the Eddington limit. 
Unfortunately, in the absence of an inclination we may not de- 
termine how close Cyg OB2 #12 is to the rotationally modified 
Eddington limit, although we note that the Eddington parame- 
ter for it is lower than Groh et al. d201 II ) found for AG Car at 
any epoch. Given that AG Car is also a known rapid rotator but 
has not been observed at, or to have exceeded, the rotationally 
modified Eddington limit we suspect the same to be true for Cyg 
OB2 #12. From this we might therefore conclude that the HD- 
limit is not the direct result of a star (b)reaching the rotation- 
ally modified Eddington limit. Indeed, Groh et al. ( 1201 II ) show 
that AG Car most closely approaches the rotationally modified 
Eddington limit during its hot, compact phases when, in contrast 
to its cool extended state, it lies beneath the HD limit (Fig. 12). 
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Appendix A: Summary of historical data for 
Galactic BHGs 

In order to understand the role of the BHG phase in the lifecy- 
cles of massive stars it is invaluable to contrast the long term 
(in)stability of such stars to other, potentially related evolution- 
ary phases such as the LBVs. In order to accomplish this we 
undertook a literature review of all known Galactic examples. 
This has allowed the construction of the comprehensive spectro- 
scopic and photometric datasets used in the quantitative model- 
ing descibed in Sect. 4. To the best of our knowledge there are 16 
known or candidate BHGs within the Galaxy, of which 8 have 
early (Bl-4) spectral types. In constructing this list we have ex- 
plicitly included 3 high luminosity stars of late-B spectral type 
which, although classified as luminosity class Iae, have promi- 
nent Balmer emission lines that appear characteristic of bona 
fide BHGs. 

A.1. The early-B HGs: Cyg OB2 #12 
A.1.1. Photometric variability 

Gottlieb & Liller ( 1978 ) provide an extensive optical photomet- 
ric history for Cyg OB2 #12, comprising, where available, an- 
nual mean B band magnitudes from ~ 1890- 1980. They report 
low level annual variability (jhb ~ 14.5 - 15.0mag.) albeit with 
no apparent secular evolution. Despite being somewhat limited, 
the compilation of multi-epoch UBVRIJHK photometry pre- 
sented in Table A. 1 is consistent with little photometric activity 
being present- 
While systematic photometric monitoring of Cyg OB2 #12 
ceased after 1980 - fortuitously at a point when spectroscopic 
observations increased in frequency (Sect. A. 1.2) - two fur- 
ther long-term photometric datasets are available. The first, ac- 
quired via the Hipparcos mission (ESA 119971 Perryman et al. 
119971 ), provides data from 1989 November to 1993 March, 
with the passband, H p , spanning 3400-9000A (van Leeuwen 
et al. |1997l l. Subsequently, the Northern Sky Variability Survey 
(NSVS; Wozniak et al. 2004 ) yielded a further year long dataset 
(1999 April - March 2000); this operated without filters, result- 
ing in a wide optical band with a response defined by that of the 
CCD, extending from 450 to 1 OOOnm. 

These data are summarised in Fig. 3 and although the non- 
standard passbands do not allow a direct comparison to previous 
data they do constrain variability over the periods in question. As 
with the results of Gottlieb & Liller (1978), these datasets indi- 
cate that low level (~0.3mag) non secular and aperiodic variabil- 
ity is present, which occurs over rather short (~days) timescales. 

While only single epoch mid-IR observations are available 
(Table A. 2), multiple epochs of radio continuum observations 
are available in the literature, and are summarised in Table A. 3. 
While the spectral indices derived from these data are consis- 
tent with thermal emission from a partially optically thick stellar 
wind, these fluxes demonstrate significant variability. As with 
the short period photometric variability, a possible interpretation 
of this behaviour is provided in Sect. 3. 



9 We note that in common with early observations of all the BHGs 
considered, interpretation is complicated by the lack of a reported date 
for some observations - a problem also present for the spectroscopic 
observations. 
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A.1.2. Spectral classification and variability 

A summary of the multiwavelength spectroscopy obtained for 
Cyg OB2 #12 over the past half century is presented in Table 
A.4. The spectral classifications given are those reported by the 
works in question; reclassification is provided only for those 
occasions where the source papers assumed an earlier result 
(this being the case for the IR studies, where Cyg OB2 #12 
was typically included as a spectral standard). Regarding 
reclassification, a luminosity class of Ia + is implictly assumed 
throughout the following discussion. 



4-6000A: Souza & Lutz (fl980l> Massey & Thompson (fl99Tb, 
Klochkova & Chentsov, d2004l > and Kiminiki et al. (f2007l > 
all provide a detailed discussion of the classification of 
Cyg OB2 #12 from 4000-5000A spectra at high resolution and 
S/N, using the Hei 4471A:Mgn 4481 A line ratio as the pri- 
mary spectral-classification criterion. Following this methodol- 
ogy, spectral types in the range B3-B8 have been inferred by 
these authors, with variability in the line ratio (e.g. Kiminiki et 
al. 120071 1 taken as evidence for changes in the spectral type of 
CygOB2#12. 

However, we caution that extreme care must be taken in 
inferring changes in the physical properties of the star from 
this diagnostic. Firstly, the Mgn 4481 A transition shows a 
dependence on stellar luminosity and so calibrations of spectral 
type versus the Hei 4471A:Mgn 4481 A line ratio derived from 
lower luminosity stars may not a priori be directly applied 
to Cyg OB2 #12. More critically and as demonstrated in 
Sect. 4.1, this ratio is highly sensitive to the properties of the 
photosphere/wind transitional zone. Consequently, changes 
in the structure of the inner wind - which we may infer from 
line profile variability in Ha (see below) - at constant stellar 
temperature will cause this line ratio to vary, and in turn lead 
to the erroneous conclusion that because the spectral type has 
varied then the temperature has also changed. 

Given the quality of their spectra - and in conjunction with 
data presented in Walborn & Fitzpatrick ( 1990) - Kiminki et al. 
(2007) also employed the weak absorption feature at ~4542A 
as a classification criterion, which they attributed to a blend 
of Hen and Fen. While an absorption feature at ~4542A is 
observed in B 1.5-4 supergiants, it is unlikely to be due to Hen 
(seen only for BO and earlier supergiants) or Fe n (B5 and later). 
Given these uncertainties, we instead choose to employ the 
Sim 4552A/Sin 4128A line ratio as a primary temperature 
diagnostic. Thus, upon consideration of the above, we suggest a 
spectral classification of B5 for the 1992 July 22 spectrum and 
B4 for the 1998 August, 2000 September and 2008 July spectra 
(Fig. A.l), noting that the marginally later classification in 1992 
is likely to be a result of the lower S/N and spectral resolution of 
these data. We thus find no evidence for changes in the spectral 
type of Cyg OB2 #12 in these data. 



6-8000A: The spectral region surrounding Ha is largely devoid 
of classification criteria (e.g. Negueruela et al. 2010). However 
(i) the strong central Ha emission peak superimposed on broad 
emission wings, (ii) lack of Hen 6527, 6683A absorption, (iii) 
absence of S iv 6668, 6701 A emission and (iv) the presence of 
strong C ii 6578, 6582A photospheric lines are uniformly present 



in spectra from 1992 onward^3 an d are all consistent with a 
highly luminous B supergiant (Fig. 1). As found by Klochkova 
& Chentsov (2004), the central peak of the Ha profile appears 
to be variable at low projected velocities (< 140kms _1 ). 



8-9000A: Despite spanning a period of over 16 years, the 
3 I band spectra are essentially identical (a representative 
spectrum from 2008 is presented in Fig. A. 2). To classify 
Cyg OB2 #12 from these data we employed the scheme de- 
scribed in Negueruela et al. (2010), supplemented by spectra of 
other early-mid BSGs and BHGs of known spectral type. While 
the BHG spectra show a similar morphological progression to 
those of normal BSGs (Negueruela et al. 120101 Fig A. 2), there 
is no indication of emission in the Paschen series (noting that 
the lower Balmer transitions are seen in emission). While such 
a finding precludes the separation of BHGs from BSGs in the 
spectroscopic window in which GAIA - the Global Astrometric 
Interferometer for Astrophysics - will operate, the lack of wind 
contamination aids in the quantitative determination of the un- 
derlying stellar parameters (Sect. 4). 

Regarding the classification of Cyg OB2 #12, the presence 
of N i absorption lines indicates a spectral type of B3 or later, 
with their comparative weakness favouring B3. This finding 
is consistent with both the strength of the Oi 8446A line, 
which clearly lies between B2.5Ia and B4Ia (see Fig. A. 2), 
and the temperature sensitive Pal5/Pall and Pal6/Pall line 
ratios. However, the differences between B3 and B4 la stars are 
small - hence we suspect that the slight discrepancies between 
classifications based on the ~4-5000A and ~8-9000A windows 
evident for data obtained in 1992 and 2008 are unlikely to be 
real. Likewise, distinguishing between B5 and B8 supergiants 
on the basis of I band spectra alone is particularly difficult 
(Negueruela et al. 2010); hence we caution against using the 
results of Sharpless ( 119571 1 and Zappala ( 1970) as evidence for 
spectroscopic variability prior to 1970. Echoing Sect. 4, changes 
in the physical properties of the photosphere/wind transitional 
zone have negligible impact on the I band spectral morphology 
- hence explaining potential discrepancies between the spectral 
type of Cyg OB2 #12 determined from the I band and the Hei 
4471A:Mgn 4481A diagnostic. 



0.98-1. lOjum: Two epochs of observations from 1990 and 1992 
are presented by Conti & Howarth ( 119991 1, who report no vari- 
ability between the spectra. Unfortunately a paucity of spectral 
features and little calibration data complicate classification 
from these data, although the presence of weak Hei 1.031yum 
suggests a classification of B8 Ia + or earlier and the absence 
of Hen 1.0124/mi provides an upper limit of B3 Ia + . Such a 
result is consistent with that derived from our spectra obtained 
on 1992 July 22. 



1.66-1.72//m: A single epoch of H band spectroscopy was 
presented by Hanson et al. (1998), with Brll and Hei 1.700yt/m 
lines in absorption and of ~equal strength and Hen 1.693/mi 
absent. The lack of Hen 1. 693pm is characteristic of a B 



The low resolution and S/N of previously published spectra are in- 
sufficient to comment on the presence, or otherwise, of photospheric 
features. 
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spectral type, with EW(Brll)>EW(Hei) occuring for B1.5 Ia + 
and later (at the resin, of the spectrum in question) and only 
becoming »1 at B8 Ia + , hence providing upper and lower 
limits to the spectral type (Hanson et al. 1998 2005 ). 



2.0-2.4yum: Hanson et al. (1996) present a K band spectrum 
from 1994, in which Bry and Hei 2.112//m are seen in absorp- 
tion. For a star of the expected luminosity of Cyg OB2 #12 and 
with the low S/N and resolution of the data in question, only 
broad constraints are possible, with the lack of Hen 2.189//m 
and the presence of Hei 2.112//m implying a classification 
between BO-9 la. 



2.35-4 .10jUm: Finally, Lenorzer et al. (2002) and Whittet et al. 
(1997) present 2.35-4. lyum spectra, which we may supplement 
with an additional spectrum from 1996 April 4. All four display 
emission in Bra, Pfa and with the higher transitions seen 
in absoprtion. No variability is apparent between any of these 
spectra; hence we do not reproduce them here. Unfortunately, 
as with the ~l-2/mi window, there is a dearth of classifica- 
tion features in this region and we may only infer a spectral 
type between B3-9 Ia + from the presence and strength of the 
Brackett and Pfund series and the lack of He i photospheric lines. 

To summarise - we find that the available data provides no 
substantive evidence for the long-term evolution of the spectral 
type of Cyg OB2 #12 over the past half century, although this 
conclusion is necessarily weaker prior to the 1990s, given the 
sparser dataset available. Nevertheless, this finding is entirely 
consistent with the long-term lightcurve. A classification of B3- 
4 Ia + appears appropiate thoughout this period, although short 
term variability is present in both wind (e.g. Ha.. .6) and photo- 
spheric lines (e.g. Kiminki et al. 120071 ) . 

Both phenomena are common in other luminous B super- 
/hypergiants (e.g. Clark et al. 120101 ) where they are assumed to 
be due to, respectively, time-dependent wind structure and pho- 
tospheric pulsations resulting in changes in photospheric tem- 
perature. Regarding the latter, changes in the Hei 4471A:Mgn 
4481 A flux ratio have led to the conclusion that the spectral 
type is variable (B3-8 Ia + ), from which changes in the photo- 
spheric temperature (16-12kK) of Cyg OB2 #12 have been in- 
ferred; such a range is fully consistent with models of the pul- 
sating B supergiant HD64760 (Kaufer et al. 120061 ) and the BHGs 
Wdl-7 and 42 (Clark et al. 120101 ). However we caution that the 
detailed non-LTE model atmosphere analysis reveals that vari- 
ation in the physical structure of the photosphere/wind transi- 
tional zone can also lead to this behaviour; hence concluding 
that the photospheric temperature has been variable based on the 
current spectroscopic dataset is clearly premature. 

A.2. The early-B HGs: C Sco 
A.2.1. Photometric variability 

Sterken (1 19771 ) and Sterken et al. (1997 ) have previously under- 
taken long-term(differential) uvby Stromgren and V BLUW pho- 
tometric campaigns between 1973-4 and 1982-95 respectively. 
They report low level (Am ~ 0.01 mag) quasi-periodic variabil- 
ity over a wide range of timescales (~ 16.5 2000 days), with 

additional stochastic variability (Am ~ 0.05mag) superimposed. 
(Aperiodic) low level (Am = 0.08mag) photometric variations 



were also present in the Hipparcos dataset between 1989-1993 
(Lefevre et al. 2009 ) and between 1979-80 (Burki et al. [19821 . 

In addition to these data Sterken et al. ( 119971 ) also present 
historical observations from the 10 -19'* Century, to which we 
may add a number of more recent data from the 20''' Century 
(Table A.l). Sterken et al. d 19971 ) suggest that Sco is a long- 
term variable, reporting my ~ 2.8 in the middle of the 18 r,! 
Century, m v -4.3 and 5.4 in 1875 and 1878 and m v ~ 4.5 - 4.8 
between 1890-1900. The latter values are broadly consistent 
with the photometric data presented in Table A.l, which span 
the 20''' Century and provide no evidence for the long term sec- 
ular variability that characterises e.g. the LBV phase. Moreover, 
these data also demonstrate that the near-IR continuum is simi- 
larly stable, albeit over a shorter, ~ 30 yr baseline. We therefore 
conclude that the sole evidence for significant (Amy > 0.3mag) 
variability is provided by historical visual estimates dating from 
before -1890, for which we are unfortunately unable to quantify 
the observational uncertainty. 

Finally, we turn to the radio data, for which Bieging et al. 
(1989 ) report possible variability. This mirrors the findings for 
Cyg OB2 #12 and we provisionally associate this behaviour with 
the same physical cause (Sect. 3). 

A.2.2. Spectroscopic classification and variability 

Given the brightness of Sco, spectroscopic observations dat- 
ing from the turn of the 19' /! Century are reported in the liter- 
ature, although a lack of accurate dates for the earlier observa- 
tions complicates their interpretation. Nevertheless, a summary 
of these observations are presented in Table A. 5; we note that 
the lower reddening to ( x Sco facilitates a more homgeneous 
dataset of blue end spectroscopy in comparison to Cyg OB2 #12. 
Fortuitously, the period between 1900-1950 is well sampled and 
hence complements the sparse photometric dataset in this period. 

This compilation reveals an absence of long-termvariability 
or secular evolution of spectral type over a ~1 10 yr interval. The 
description of the main features of the optical spectra over this 
whole period are remarkably consistent; the spectra presented 
in Fig. A.l for the period 1994-2009 demonstrating this stabil- 
ity. Nevertheless, numerous authors report line profile variabil- 
ity in the wind dominated P Cygni profile of Ho- and with 
the higher H i transitions repeatedly varying between emission 
and absorption over at least a ~50 yr interval (cf. Table A. 5 and 
refs. therein). Rivinius et al. (1997) report the rapid variability 
of both photospheric and wind lines between 1992-5; we con- 
sider it likely that the variability observed in the Ho- P Cygni 
line between 1998-2009 (Fig. 1) reflects a continuation of this 
behaviour. 

A.3. The early-B HGs: HD 190603 

Unfortunately, photometric datasets for HD 190603 are some- 
what sparse in comparison to the previous two stars, with no 
photometric observations reported over the last decade and a 
>16 yr gap in the 1980s-90s (Table A.l). Nevertheless, we find 
no evidence for secular variability over the ~46 yr period from 
1952-1998. Conversely, rapid photometric variability has been 
reported on several occassions by both Percy & Welch (1983) 
and Koen & Eyer d2002l ). 

Likewise comparatively few spectral observations have been 
reported, a problem exacerbated by the lack of reduced spectra 
being presented. Consequently we are forced to simply present 
the spectral classifications reported for HD 190603 in Table A. 6. 
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These data indicate a corresponding lack of spectroscopic evo- 
lution over a 50 year period. Indeed, the description of the spec- 
trum by both Beals d 19461 1 and Merrill & Birwell d 19431 ) as, re- 
spectively, that of a P Cygni supergiant and a Be star, suggest a 
similar morphology prior to 1933. The latter authors further de- 
scribed variable emission lines; similar findings being reported 
by Rosendhal (119731 1 and Rivinius et al. (11997I ). 

A.4. The early-B HGs: HD 80077, HD 169454 and BP Cru 

We next turn to the remaining early-B HGs. We refrain from tab- 
ulating their sparse photometric datasets but note that there is no 
evidence for secular variability in any of the three stars; for in- 
stance comparison of the optical (Hiltner [1956l Kilkenny et al. 
[19931 and near-IR (Whittet et al. [19971 Skrutskie et al. 120061 
data for HD 169454 reveal constancy over ~35 and ~26 yr inter- 
vals respectively. However, following the previous discussions, 
rapid, low amplitude photometric variability appears ubiquitous 
(van Leeuwen et al. 119981 Sterken 119771 and Hammerschlag- 
Hensberge et al. 1979 respectively). 

As with HD 190603, Table A. 6 summarises their spectral 
types as reported in the literature. No evidence for the long-term 
(secular) evolution of spectral type over timescales in excess of 
30 yrs was found, with the description of HD 169454 as a B 
star with Hi P Cygni emission lines (Merrill & Burwell 1943) 
suggesting spectral stability for nearly a century. In contrast line 
profile variability (LPV) on a timescale of days-weeks is present 
in all 3 star s (Kno echel & Moffat [T982l Rivinius et al. [19971 and 
Kaper et al.[2006 respectively). 

A.5. The early-BHGs/WNVL stars Wd1-5 and 13 

The remaining early-B HGs are found within the massive young 
cluster Wdl, and on the basis of a restricted wavelength range 
(5800-8900A) were classified as borderline BHG/very late WN 
stars, forming an evolutionary sequence from Wdl -5 through - 
13 to the WN9h star Wdl-44 (Clark et al. l2008l . Unfortuntely, 
these stars have only been observed over the past -decade, but 
no evidence for long term variability has been found (Clark etal. 
2010). However, Wdl-13 is a confirmed -9.27 day massive bi- 
nary, while pronounced LPV in Wdl-44 also argues for such 
an identification. Consequently, we suspect that all three stars to 
have formed as the result of close binary evolution. 

A.6. The late-B HGs 

Finally, we examine the eight BHGs with spectral types of B5 
and later that have been identified within the Galaxy. Three - 
W7, W33 (both B5 Ia + ) and W42a (B9 Ia + ) - are located within 
the massive young cluster Wd 1 and are discussed in detail in 
Clark et al. (120101 ). Unfortunately, the long-termspectroscopic 
and photometric datasets for these stars are less complete than 
those of the early-B HGs described above, although they are suf- 
ficient to confirm the presence of rapid LPV. This behaviour, as 
well as short term low amplitude photometric pulsations is also 
present in the other examples; HD 160529 (A9-B8 Ia + ; Stahl et 
al. 120031 , HD 168607 (B9 Ia + ; Chentsov et al.l2003l Sterken et 
al. [19991 . HD 168625 (B8 Ia + ; Sterken et al. [19991 Chentsov et 
al. 120031 , HD 183143 (B7 Iae; Adelman et al. l2000l Chentsov et 
al.HUO l and HP 199478 (B8 Iae Percy et al. l2008l Markova & 
Valchev 2000), suggesting that the a Cygni instabilities are ubiq- 
uitous across the complete temperature and luminosity range 
spanned by BHGs. 



This phenomenon is also found to extend to cooler temper- 
atures, having been identified in a number of early-A (AO-2.5 
Ia + /Iae) stars with similar spectral morphologies to the BHGs 
- e.g. HD 92207 (AO Iae; Sterken [19771 Kauferet al. [1997l l. 
HD 2 23385 (A2.5 Ia + ; Adelman & Albay rak [19971 Chentsov 
120031 and AS 314 (AO Ia + ; Miroshnichenko l2000l . 

However, unlike the early-B HGs, HD 160529 and 168607 
demonstrate characteristic LBV photometric modulation, with 
the former also exhibiting correlated spectral type variability, 
while both HD 168625 and HD92207 show evidence for a com- 
plex, dusty circumstellar environment (Roberto & Herbst 1998; 
Clarke et al. 2005 ) implying recent enhanced mass loss, possibly 
associated with an LBV phase. 
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Fig. A.l. Long-term evolution of Cyg OB2 #12 (1992-2008; top panel) and £ 1 Sco (1998-2009; bottom panel) in the -4000-4800A 
spectral region, encompassing the Sim 4552A/Sin 4128A temperature diagnostic. 
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Fig. A.2. Comparison of the I band spectra of Cyg OB2 #12 and Sco with other galactic BHGs. Note the systematic weakening 
of He i and strengthening of N i transitions with decreasing temperature, which permit spectral classification in whis wavelength 
range. 
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Table A.l. Summary of historical photometric magnitudes used to construct the optical-near IR SEDs, ordered by time. 



Date 








Waveband 








Reference 




U 


B 


V 


R I 


J 


H 


K 














Cyg OB2 #12 












<1954 


16.67±0.1 


14.66±0.04 


11.47±0.04 


- 


- 


- 


- 


Sharpless (19571 


<1962 


17.21 


14.71 


11.49 


8.36 6.20 


4.55 


- 


2.81 


Johnson & Borgman 


1963) 


<1964 


17.20±0.01 


14.70±0.01 


11.48±0.01 


8.26±0.01 5.94±0.01 


4.32±0.01 


_ 


2.66±0.1 


Johnson (1965 1 


<1967 


17.15 


14.70 


11.48 


8.26 5.95 


4.42 


_ 


2.70 


Wisniewski et al. (119671) 


<1968 


_ 


_ 


_ 


8.26 5.94 


4.32 


_ 


2.66 


Johnson (1968) 


<1971 


_ 




_ 




- 




2.60 


Stein & Gillett ( 1971 1 


<1973 


- 




- 




_ 


3.34±0.06 


2.77±0.05 


Chalduetal. (1973) 


1972 Autumn 


- 


- 


- 


- 


4.59±0.05 


3.37±0.05 


2.82±0.05 


Voelcker dl975l) 


1975-6 










A TO i A AO 

4.38±0.02 


3.28±0.02 


2. /z±0.0z 


Harris et al. (11978b 


1980 Aug 26 








_ 






2.80±0.01 


Abbott etal. 0981) 


1981 July 


_ 


- 




- 


_ 


3.33±0.01 


2.72±0.01 


Leitherer et al. (11982b 


1989 May 26-30 


- 


- 


- 


- 


4.39±0.04 


3.34±0.04 


2.75±0.04 


Torres-Dodgen et al. (1991 


1990 June 6 


17.18+0.01 


14.81±0.01 


11.46±0.01 


- 


- 


- 


- 


Massey & Thompson ( 


1991) 


1998 June 21 


- 


- 


- 


- 


4.7±0.3 


3.5±0.3 


2.70±0.4 


Skrutskie et al.( 2006 ) 










C Sco 












<1924 


_ 


- 


4.88 




_ 


- 


- 


HD catalogue 


1949- 1952 1 






4.88 










Schilt& Jackson (|1952|( 


<1954 2 


_ 


_ 


4.9 


_ 


_ 


_ 


_ 


Bidelman d 1954b 


<1958 


_ 


_ 


4.86 


_ 


_ 


_ 


_ 


Code &Houck (1958) 


<1958 


_ 


5.03 


4.3 




_ 






Evans et al. (1959) 


1958 March- April 


4.54 


5.20 


4.82 


_ 


_ 


_ 


_ 


Westerlund (fl958l> 


<1961 




_ 


4.80 










Buscombe (119621) 


1963 


_ 


_ 


4.74 










Andrews et al. ( 1968) 


1963 June 


4.57 


5.19 


4.72 


- 


- 


- 


- 


Feinstein (119681) 


1964 March 


4.64 


5.20 


4.72 


- 


- 


- 


- 


Feinstein dl968b 


1967-68 


- 


- 


4.71 


- 


- 


- 


- 


Schild et al. (|1969|( 


1968 July 


- 


- 


- 


- 


- 


3.30 


3.11 


Schild etal. (1971) 


<1973 


- 


- 


4.73< V <4.83 


- 


- 


- 


. 


Jaschek& Jaschek dl973|( 


1974 April 


- 




- 


- 


3.47±0.03 


3.27±0.03 


3.11+0.03 


Whittet etal. (1997) 


1963 to 1976 3 






4.73 










Feinstein & Marraco 


1979) 


1979 Apr. to 1980 Sept 4 


- 


- 


4.66< V <4.78 


- 


- 


- 




Burki et al. 


1982) 


i no 1 a . i~\ 

1981 August 9 














i 1 A i ft r\c 

3.14±0.05 


Abbott et al. dl98lb 


1974-82 


4.76 


5.29 


4.77 


4.31 3.85 




_ 




Fernie 


1983) 


1982 August 


_ 




_ 




3.58±0.03 


3.32±0.03 


3.18±0.03 


Leitherer & Wolf ( 


1984) 


<1984 










3.53±0.02 


3.31±0.02 


3.13±0.03 


Lopez & Walsh 


1984) 


<2003 


4.61 


5.18 


4.78 


_ 


_ 


_ 


_ 


Reed et al. 


2003) 


1999 May 22 


„ 


_ 


„ 


_ 


3.6±0.3 


3.3±0.2 


3.3+0.2 


Skrutskie et al.( 2006) 










HD190603 












1952-4 


5.72 


6.18 


5.65 










Hiltner ( 


1956b 


<1958 


5.59 


6.13 


5.56 










Golay ( 


1958) 


<1961 




6.14 


5.59 










Bouigue et al. ( 


1961) 


1962-70 


5.66 


6.14 


5.62 










Crawford et al. 


1971) 


<1964 


5.72±0.01 


6.18+0.01 


5.65+0.01 


5.12±0.01 4.72±0.01 


4.50±0.01 




4.11+0.1 


Johnson 


79651) 



Table A.l. continued. 



Date 


U 


B 


V 


Waveband 
R I J 


H 


K 


Reference 


<1966 


5.73 


6.19 


5.65 




- 


- 


Johnson et al. (1996) 


<1968 


5.69 


6.15 


5.60 


- 


- 


- 


Blanco et al. dl968l> 


<1968 


- 


6.14 


5.62 




- 


- 


Lesh(fT968l 


1968-9 


5.81 


6.25 


5.70 


- 


- 


- 


Burnichon & Gamier ( 19761 


1972 Autumn 


- 


- 


- 


4.50±0.05 


4.14±0.05 


4.11±0.05 


Voelcker dl975l> 


<1974 


5.72 


6.18 


5.64 


- 


- 


- 


Nicolet(1978) 


1974 Sept. 








4.43±0.03 


4.24±0.03 


4.10±0.03 


van Breda & Whittet ( 1977 1 


1974-82 


5.64 


6.13 


5.56 


5.01 4.62 






Fernie (19831 


1975 Oct. 












4.1+0.2 


Sneden et al. (19781 


1976-78 




6.18±0.02 


5.63±0.02 


5.08±0.02 4.71±0.02 






Moffett & Barnes (fl979l> 


1979-82 




6.18 


5.64 








Percy & Welch (fl983ll 


1998 May 11 








4.42+0.26 


4.26±0.20 


4.04±0.03 


Skrutskie et al. (12006b 



Note that where available errors and exact dates of observations are presented. For brevity we have not replicated the B band lightcurve of Gottlieb & Liller ( 19781 for Cyg OB2#12 here. 
Observations between 1945-1951 by Cousins d!952t suggest that ( 1 Sco was variable with low amplitude (Am = 0.05mag) between 1945-51, although details of the filter set used are not described. 
Sterken dl977t report 35 epochs of differential ubvy photometry between 1973-1974 for Sco, with a further 399 epochs between 1982-95 presented by Sterken et al. (1997); throughout these 
periods minimal variability was observed (A ubvy <0. lmag.). Finally Percy & Welch ( 1983 1 report 52 epochs of BV band photometry between 1979 May 30 and 1982 July 1 indicating t to be a 
short period variable. '3 observations in this period; no details of band pass 2 Phot. magnitude possibly sourced from HD catalogue? 3 3 observations over this period. 4 43 observations over this 
period. 
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Table A.2. Summary of mid- to far-IR fluxes used to construct 
the SEDs 



Wavelength 


Magnitude 


Reference 


Qim) 








Cyg OB2 #12 


2.30 


2.72 ±0.01 


Leitherer et al. (1982) 


2.345 


2.635±0.02 


Harris et al. {1978 > 


3.013 


2.405±0.02 


Harris et al. (11978b 


3.108 


2.375±0.02 


Harris et al. (1197811 


3.420 


2.281 ±0.02 


Harris et al. (11978b 


3.45 


2.37±0.05 


Abbott et al. (T98TT) 


3.50 


2.28±0.05 


Chalduet al. (1973) 


3.57 


2.28±0.02 


Leitherer et al. dl982b 


3.57 


2.170±0.05 


Torres-Dodgen et al. ( 1991 ) 


3.58 


2.217±0.02 


Harris et al. (1978) 


3.820 


2.178±0.02 


Harris et al. (1978) 


4.29 


2.057±0.118 


MSX (Eganetal. 120011) 


4.80 


2.37±0.10 


Abbott etal. (1981) 


4.97 


2.06±0.02 


Leitherer et al. (11982b 


8.28 


1.946±0.044 


MSX (Eganetal.l200TT> 


10.2 


2.03±0.10 


Abbott etal. (1981) 


10.9 


1.95±0.07 


Leitherer et al. (1982) 


12.0 


1.642±0.00 


IRAS PSC (1985) 


12.13 


1.900±0.056 


MSX (Eganetal.l200lb 


14.65 


1.713±0.067 


MSX (Eganetal. 120011) 


20.0 


1.54±0.10 


Abbott etal. (1981) 


21.34 


1.575±0.079 


MSX (Eganetal.p200n> 


25.0 


1.204 


IRAS PSC (1985) 



C Sco 



3.45(L) 2.98±0.03 Leitherer & Wolf ( 1984 ) 

4.7(M) 2.90±0.03 Leitherer & Wolf (1984) 

8.28 2.624 MSX (Egan et al. l200l1 ) 

12.13 3.416 MSX (Egan et al. 120011 ) 

14.65 3.560 MSX (Egan et al. [200TT > 



3.6 

4.9 

8.28 

9.0 

12.13 

14.65 



3.96: 

3.89±0.20 
3.68±0.04 
3.12+0.01 
3.47±0.08 
3.46±0.09 



HD190603 

0.18 



Sneden et al. (19781 ) 
Snedenet al. (1978) 
MSX (Egan etal. 120011 ) 
AKARI 
MSX (Egan etal. 120011 ) 
MSX (Egan etal. 120011) 
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Table A.3. Summary of mm & radio fluxes in mJy used to construct the SEDs, ordered by time. 



Date 


u. ijcm 


O 1 nm 

v. 1 cm 


1.3cm 


2cm 


Waveband 

3.5cm 3.6cm 


6. leni 


6. 2cm 


70 S,-ni 

zu. jem 


Reference 












Cyg \JdZ tfil 










lyov) Marcn ZZ 






Zl±o 












WenuKer &l Altennoti ( lysu ) 


1Q0fl \/Tn.\r i 01 

ivou May zz+zj 














o.Z+U.j 




ADDott et ai. \ iVoi ) 


IVol UCt. ID 








O.U+Z.U 






1 a _i_n 7 
0.4+U.Z 




Dieging et ai. I lyoy ) 


iVoz Marcn z / 


















wmte & KecKer i rvoJ i 


1982 Aug. 26 


_ 


_, 


_ 




_ 


_ 


2.2+0.2 


2.0±0.2 


Bieging et al. (1989) 


1983 Aug./Oct. 


- 


- 


- 


- 




4.0+0.3 


- 


2.1±0.3 


Becker & White ( 1985 ) 


I77J May 1 










Of^-i-0 07 
O.UD+U.U / 








waiaron ei ai. \ vyyo i 


1 QQ4. Spnt 

l77 i T OCpi. 








7 7+n ^ 


4. 74+0 14 
*+. / tiu. it 




4 00+0 70 




oluucii cl ai. niyyot 


tW4 UCl. 








1 7 n_i_n 7 

tZ.U+U.Z 


/.4U+U.U0 




/i nn_i_n 7n 
4.UU+U.ZU 




ocuuen et ai. \ two i 


177J /\pm z / 




77 Q_i_n 

ZZ. V+U.O 




1 1 3-1-0 1 


/.io+u.U4 


J.04+U. 1Z 






i^onireras ei ai. ( ryyo i 


1 QOQ Iiino 7Q 

ivy? June Zo 




77 Q_i_n a 
zz. y+u.o 




n_i_ 1 ^ 


^ Q_i_n 1 


a 7_i_n 1 

4.Z+U. 1 






contreras et ai. ( zuu4 1 












<r' sco 










IVol May / 












l.Z+U.J 






riieging et al. ( Ly<$y ) 


1V54 Marcn y 












7 n_i_n 7 
Z.U+U.Z 






Dlarriurt at- nl d 10001 

riieging et al. fliyoyi 


1984 April 3 








4.3+0.1 




1.9+0.2 






Bie^inp et al dl 9891 


1990 Sept. 12-14 


23.0+2.4 
















Leitherer & Robert d!991t 


<2006 










2.4 








Benaelia et al. (120071 












HD 190603 










1979 July 13 












<0.5 






Abbott etal. (|198U 


1994 Oct. 








0.7+0.2 


0.7±0.1 




0.58+0.06 




Scuderi et al. d 1998 \ 



Note that for brevity the wavelengths are presented to the nearest 0.1cm, with exact values given in the papers in question. 
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Table A.4. Summmary of historical spectral observations of Cyg OB2 #12 



Date 


Wavelength 


Spectral Type 


Comments 


Reference 


<1954 


4100-6600A 


<A0I 


Ha in emission 


Morgan et al. ( 1954) 


<1957 


8-9000A 


B51a 


- 


Sharpless 4.1957} 


<1970 


near IR 


B8Ia + 


- 


Zappala (|l970|l 


<1971 


? 


- 


Ha in emission 


Bromage 4 197 111 


<1973 


5-8700A+ 


- 


Ha in emission 


Chalduetal. (|1973|| 




8-11000A 








1977 July 


3840-4480A 


B8 Ia + 


Variable Ho- 


Souza&Lutz(|l980|l 




+Ha 




rn emission 




1980 Autumn 


Ha 


- 


Ha in emission 


Leithereret al. 41982b 


<1981 


4000-6600A 


- 


Ha+HB in emission 


Hutchings (|l98l|l 


1989-90 summer + 


encompasses Ha 


B5 la 


Ha+HB in emission, 


Massev & Thompson 4 199 111 


'few years earlier' 


+ -4089-4686A 


B5 la 


Hy+H6 infilled 




1990 August 


0.98-1. 10/zm 


>B3 Ia,<B8 la 


- 


Conti&Howarth(|l999|l 


1992 July 


0.98- 1.10pm 


>B3 la, <B8 la 




Conti&Howarth(|l999|l 


1992 July 17 


4036-4836A 


B4-5 la 


11 t H in absorption 


This work 




5765-9586A 


B3 la 


Ho- in emission 




1994 Sept. 


2-2. 2pm 


>B0.5 la, <B9 la 




Hanson et al. (|1996|) 


1995 July 17 


6366-6772A 


- 


II in emission 


This work 


1995 Dec. 23 


2.35-8pm 


>B3 la, <B9 la 


Bra-, Pfa+B in emission 


Whittetetal. (|1997|| 


1996 April 4 


2.35-8pm 


>B3 la, <B9 la 


Bra, Pfa+B in emission 


This work 


1 QQ^ f V't 1 7 

iyyo wet. i / 


z.jj-o//m 


>Dj la, <Dy la 


Bra, Pfa+p 1 in emission 


wnittet et ai. u lyy / p 


1 997 Tune 


1 <S(S-1 7?/;m 

1 . VJlJ 1 . / /_ fill 


^ 1_> 1 . .. ' la, NUO 1 tl 




Hanson et al it 1 998h 

ii.aiistjii c-l ai. \ 1770 I 


1998 Apr-May 


2.35-8pm 


>B3 la, <B9 la 


Bra, Pfa+B in emission 


Lenorzer et al. (2002) 


1998 Aug. 09 


3940-5737A 


B4 la 


Hp 1 in emission 


This work 




6366-6772A 




Ha in emission 




iyyy July 1U 


jjUU- / /UUA 


rS la 


Ha in emission 


Kiminki (2010, priv. comm.) 


ZUUU July 1U 


jjUU- / /UUA 


rS la 


Ha in emission 


Kiminki (2010, priv. comm.) 


oooo c nt 1 c 
zUUU oept. 1 o 


TAOO ^OOOA 
:>oUUOZUUA 


DQ Too 

rJj lae 


np inniieci 


xviminKi ei ai. (zuu/ ) 


1 AAA 4 11 

2000 Sept. 25 


/I AAA ylTCA A 

4000-4750A 


B4 la 


Hy+H(5 in absorption 


This work 




6340-6740A 




Ha in emission 




2001 June 21 


4542-7939A 


B5 la 


B aimer line emission 


Klochkova & Chenstov fl2004b 


9001 Anmict 9zL 
ZUU1 rVUgUSL Zf 


1SOO ASOOA 


DA I.,- 
DO IdC 


Hy+H(5 in absorption 


JVlIIllIlKJ eL al. \ZAj\jI ) 


2001 Sept. 01 


3800-4500A 


B8 lae 


Hy+H<S in absorption 


Kiminki et al. \2QQl\ 


2003 April 12 


5273-6764A 


B5 la 


Ha emission 


Klochkova & Chenstov 42004|l 


2007 Aug. 22 


7600-8900A 


B3 la 




This work 


2007 Aug. 30 


5600-6800A 


B la 


Ha in emission 


Kiminki (2010, priv. comm.) 


2008 July 22 


4000-5270A 


B4 la 


1 1/3 in emission 


This work 




6450-7150A 




Hy+H(5 in absorption 








Ha in emission 






8350-8900A 


B3 la 







Note that a parameter is listed in italics if not explicitly given (wavelength of observation) or the data from which it is derived was not presented. 
Parameters given in bold are from this work. Ha is variable between the two epochs of observations presented by Klochkova & Chenstov ( 2004 1 
as Hp" and y are between the 2000-8 spectra presented here, while Conti & Howarth ( 1999 1 report no variability between 1990-2 in the ~ ljim 
region. Selected spectra between 1992-2008 are plotted in Fig. A.l. 
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Table A.5. Summary of historical spectral observations of f 1 Sco 



Date 


Wavelength 


Spectral 
Type 


Comments 


Reference 


<1896 


- 




HB in emission 


Pickering (18961 


1891-1899 


4000-4900A 


Blp 


P Cygni profiles in HB, y 
all other lines in absorption 
t Hp i 4009 On 4149 fn 4767) 


Cannon jl90l|t 


1924 


4100-4900A 


B(2)eq 


P Cygni profiles in HB, y 

1 1U III UUlIIJi LtllAJIl 


Merrill et al. jl925|t 


1903-1929 




Blela 




Bidelman jl954|) 


1QOA 1 QQA 
LyZO- LyjVJ 


jyju-'i-yUUA. 


tsipe 




Kimmer uiroup 


30/07/34 






P Cygni profiles in Ha and D3 He I 


Merril & Burwell (1943b 


<1938 






Listed as variable 


Payne-Gaposchkin & 
Gaposchkin (1938) 


1954-58 


3800-6600A 


B1.5 Ia + 


No appreciable difference 
from Cannon 1190U 


Code &Houck( 1958) 


26/04/58 to 03/09/60 1 


3900-4500A 


Blela 




Buscombe et al. (1962|) 


22/08/61 


6200-6600A 




Ha in emission 


Jaschek et al. 11 9641) 


7/08/66 


3800-4700A 


Bl Iae 


HB in emission, Hy, S 
Or He I 4471 in absn. 


Buscombe & Kennedy 
{19681) 


6/07/60 

9/08/60 

10/08/60 

28/06/66 

27/07/66 


3100-6750A 


Bl la 


P Cygni profiles in 
Ho-, B&He M5876, 6678 
Remaining lines in absorption 
Variability in Ha profile, but 
no global spectral evolution 


Hutchings ||1968|) 


<1969 


blue end 


B1.5Ia 




Schild etal. (1969) 


21/04/65, 20/05/67 
18/03/68, 14/08/70 


3500-4900A 


Bl Ia + 


P Cygni profiles in HB, y 
Variable He i, N n & n phot, lines 


Jaschek & Jaschek 
(19731 


6/08/68 


5800-6700A 




P Cygni profiles in 
Ha, He lU5876, 6678 


Rosendhal jl973|) 


16/05/73 to 10/07/73 


blue end 




No apparent variability in spec, type 


Hutchines et al. 11976b 


28/04/74 to 1/05/74 
28/04/75 


3800-4400A 


B1.5Ia + 


No variability above 20% in 
line strength reported 


Walborn jl975|t 


1974 Oct. & 1975 Mar. 


~Ha 




Variation in Ho- P Cygni profile 


Dachsetal. jl977|) 


26/06/72 to 14/08/75 2 


3400-6700A 


Blla 


Variable P Cygni profiles in 
Ha,j8& He 1^5876, 6678 
No secular evolution in spec, type 


Sterken & Wolf 
jl978|) 




5500-6800A 


d i . j la 


VdilaUlC r \^.yglll piUlllCa 111 

Hff,/3&He 1/16678 


IjUpCZ, tx, VValMl 


1 QQC4 

Lyyj 

1990 to 1994 3 


4000-6740A 


R 1 ^ To + 

dI.j la 


Variable P Cygni profiles in 
Ha,B, Heii6678 &Fem 
Hy, S in absorption 
No changes in He nSim ratio 
& hence spec, type 


Diiiinini' nt ol lMOO'7ll 

Kivinius et al. ( iyy / ) 


lyyH Oct. 7 


4UUU-895UA 


B1.5 Ia + 


■ ■ i , 11 | //wn n / ■ . . . " 1 ■ 1 

Ha, p, He i/l 6678 Y Cygni, higher 
Balmer series in absorption 


This work 


1999 June-July 


/I AAA OACA A 

4UUU-895UA 


B1.5 Ia + 


■ ■ II 11 I / /' -T 11 1 1 1 • . . .i ' 1 ' 1 

Ha, p, He lA 6678 r Cygni, higher 
Balmer series in absorption 


This work 


2002 March-May 6 


5810-7205A 




Variable Ha P Cygni 


Morel et al. ( 2004 ) 


2003 May- June 


4000-7000A 


B1.5 Ia + 




Crowther et al. (2006|t 


2005 April 24 


3800-6800A 


B1.5 Ia + 


Ha,B & He IJ6678 P Cygni, higher 
Balmer series in absorption 


This work 


2006 Feb. 16 


3933-7985A 


B1.5 Ia + 


Hafi & He 1/16678 P Cygni, higher 
Balmer series in absorption 


This work 
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Table A.5. continued. 



Date 


Wavelength 


Spectral 
Type 


Comments 


Reference 


2009 March 06 


3060-5600A 


B1.5 Ia + 


US P Cygni, higher Balmer 
series in absorption 


This work 



Note that a parameter is listed in italics if not explicitly given (wavelength of observation) or the data from which it is derived was not presented. 
Selected spectra between 1992-2009 are plotted in Fig. A.l. '6 observations in this period, 2 44 observations in this period, 3 233 observations in 
this period, 4 57 observations in this period, 5 14 observations in this period. 6 6 observations in this period. 
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Table A.6. Summary of the dates of historic spectral classifica- 
tions of the Galactic BHGs. 



Publication 


Observation 




Date 


HD 80077 (B2.5 Ia + ) 




Morgan et al. (119551 


1942-55° 


Buscombe & Kennedy (119681 


1963-67* 


Nordstrom (119751 


<1975" 


Moffat & Fitzgerald d 1977k 


<1977" 


Knoechel & Moffat (|1982) 


1977 a,6 


Negueruela et al (in prep.) 


2008 


HD 169454 (Bl Ia + ) 




Morgan et al. (119551 


1942-55" 


Hiltnerd 19561 


<1956 


Code&Houckdl9581 


<1958 


Botto & Hack d 19621 


<1962 


Hutchings J19701 


1968 


Wolf & Stahl JT9851 


1972-4 


Walborn (fT980t 


1973-5 


Hutchings d 19761 


1973-5 


Rivinius et al. d 19971 


1992-5* 


Hanson et al. dl9961 


1994 1 ' 


Hanson et al. dl998k 


1997'' 


Groh et al. <|2007) 


2001-4" 


HD 190603 (B1.5 la ) 




Morgan et al. dl9551 


1942-55 


Ahmad <|1952|> 


<1952 


Hiltner (fT956t 


<1956 


Slettebak (fT956]> 


<1956 


Lesh {19681 


<1968 


Hutchings fl970l 


1968 


Walborn dl97U 


1969 


Hutchings J19761 


1973-5 


Bisiacchi et al. ( 1976) 


1975-6 


Andrillatet al. {1995) 


1988-93 


Rivinius et al. 997 1 


1990-1* 


Lennon et al. dl9921 


1990 


Hanson et al. dl996k 


1994^ 


Blumet al. 41997b 


1996^ 


Hanson et al. (119981) 


1998 c 


Crowther et al. (2006 ) 


1992-2003 


Markova & Puis d2008|) 


<2007 


BP Cru (Bl Ia + ) 




Vldai ( 19 IS I 


19/ J 


Bord(fT979k 


1974 


Hammerschlag-Hensberee et al. d 19791 


1975 




1977 




1978 


Parkes et al. \19&0\ 


1977-78"' 


Kaper et al. dl995b 


1984 


Kaper et al. d20061 


1996* 




2002* 



All classifications made in the optical(4-6000A) band unless otherwise 
noted - "marginally earlier spectral type of B2 Iae given. *multiple spec- 
tra in this period. Observations within the near-IR (~l-2.2yum) window, 
"marginally later spectral type of B 1.5 or B2 Iae given. 



Appendix B: Spectropolarimetry 

Considering a simple 'core-halo' wind model for heuristic pur- 
poses, electron scattering of photospheric radiation in an asym- 
metric outflow will generate a grey intrinsic linear polarization, 
while emission lines formed in the wind will see a smaller scat- 



tering optical depth, and so are expected to be less polarized 
(e.g., McLean 1 19791 ). This will result in depolarization through 
the line, although the addition of an interstellar-polarization vec- 
tor means that this 'line effect' often manifests in other ways, 
including an observed increase in degree of polarization, P. 

Previously unpublished spectropolarimetric observations of 
Cyg OB2 #12 and Sco were obtained as part of the in- 
vestigations reported by Harries, Howarth & Evans (2002). 
Cyg OB2 #12 was observed on 1995 July 17 (using the WHT 
with ISIS spectrograph; R 3300), and \ x Sco on 1997 Jun 5 
(AAT with RGO spectrograph; R 5000), with the data acqui- 
sition and reduction essentially as in the manner described by 
Harries & Howarth ( 1 19961 ). Results are displayed in Fig. B.l. 

For Cyg OB2 #12 there is no compelling evidence of a 
change in polarization through the Ha emission line (although 
there is a hint of a possible increase in P), and hence no strong in- 
dication of large-scale asymmetry in the electron-scattering en- 
velope. Similarly, there is no strong evidence for temporal po- 
larimetric variability; although the spread in published /?-band 
photopolarimetric measurements is somewhat larger than their 
quoted formal errors, the results presented by Schmidt et al. 
( [19921 P = 7.893+0.037%, 6 = 116.23 + 0.14°) Whittet et 
al. ([19921 7.97+0.05%, 1 17°), and Kobulnicky, Molnar & Jones 
([T9941 8.35+0.21%, 117.8°) are all broadly consistent with the 
present results, as are earlier O-band measurements reported by 
Kruszewski (119711 ) and by McMillan & Tapia ( 119771 ). We note, 
however, an anomalous result reported by Schulz & Lenzen 
(TT9831 P = 10.05+0.12%, 9 = 125°); their observations of 
two further stars agree well with measurements presented by 
Schmidt et al. (fT992l ). 

In spite of its brightness, ( l Sco has been much less exten- 
sively observed, with the only available point of direct compar- 
ison being the narrow-band (AA = 8.5A) Ha polarimetric mea- 
surement of P = 2.39 + 0.03%, 6 = 64.6° reported by McLean 
& Clarke (1979)0 This is in agreement with our results; how- 
ever, our spectropolarimetry clearly reveals a significant, ~ 1.5° 
position-angle rotation through the emission line. 

Physically, the line effect may be associated with large-scale, 
axisymmetric structures (such as might result from rapid rota- 
tion), or irregular 'clumps'. For a star with parameters similar 
to those of Sco, summarized in Table 3, the critical equato- 
rial rotation velocitjQ is v e 210 km s _1 . With a measured 
v e sin i of ~45 km s _1 (Table 3), the intrinsic rotation is probably 
substantially subcritical (v e /v cr ; t < 0.69 with 95% confidence). 
A rotationally-induced axisymmetric departure from spherical 
symmetry (i.e. a 'disc' or 'polar' wind) therefore seems a pri- 
ori improbable as the origin of the line effect in this star. A 
more likely cause of the observed change in polarization is the 
presence of transient large-scale inhomogenities, or 'clumps', 
as has been proposed for P Cygni and AG Car on the basis of 
spectropolarimetric observations by e.g. Nordsieck et al. (2001) 
and Davies et al. (|2005); indeed CMFGEN model atmosphere 
analysis of both these stars confirms the presence of the signif- 
icant wind clumping (Najarro et al. 120011 Groh et al. [2009a). 
Therefore, by direct analogy, an additional observational test 
of this hypothesis for f 1 Sco would be the detection of time- 
dependent variations in polarization (under the assumption that 
there is no prefered geometry for the clumping), noting that 



11 Matthewson & Ford ( 1970 ) report a blue-band measurement of P = 
2.29%, e = 67.3°.) 

12 The velocity at which the outward centrifugal force equals the in- 
ward gravitational force. 
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Fig. B.l. Linear spectropolarimetry of Cyg OB2 #12 and f 1 Sco. 
Rectified spectra are shown at full resolution; polarimetry is 
binned to an (almost) constant error of 0.05% in degree of polar- 
ization. 



spectroscopic observations are already strongly indicative of the 
presence of transient wind structure (Appendix A). 

Appendix C: The age of NGC 6231 and Sco OB1 



Given its proximity (~1.64kpc; Sana et al. 2008), NGC 6231, its 
host association Sco OBI, and the massive stars located within 
both have been the subject of numerous multiwavelength studies. 
Here we review these data in order to constrain the evolutionary 
history of both cluster and association to better understand the 
nature of Sco. Photometric studies by Baume et al. £1999) 
suggest a significant age spread within NGC6231, with star for- 
mation apparently commencing ~ 10 Myr ago, culminating in 
the formation of the massive stellar cohort ~ 3 - 5 Myr ago. The 
latter age is consistent with the findings of Sana et al. (12006a, 
120071 ). who report a similar value of 2 - 4 Myr from an analy- 
sis of the photometric data of (X-ray selected) pre-MS stars. An 
analagous study of the wider Sco OB 1 association by Perry et 
al. 0991) supports a non coeval star formation history for this 
region, yielding an age of ~ 8^3 5 Myr. 

Catalogues of the OB stellar conten t of bo th NGC 6231 and 
Sco OBI are provided by Sana et al. (2006b) and Ankay et al. 
(2001) respectively, and allow for the individual placement of 



stars on the HR diagram presented in Fig. 13. For those stars in 
Sana et al. we adopted the luminosities given by these authors 
and utilised the spectral type/temperature relation of Martins 
et al. d20051 l. The temperatures and luminosities of the stars 
in Ankay et al. were again determined via the callibrations of 
Martins et al., with reddening for individual stars calculated via 
the intrinsic colours of Martins & Plez (2006). Finally, given 
their rather evolved nature, we adopted the results of the tailored 
non-LTE analysis of Crowther & Evans ( 120091 for HD 151804 
and 152408. 

Comparison of these data to non-rotating and rotating 
Geneva isochrones clearly indicate that the region as a whole 
appears non coeval, at first glance being consistent with a spread 
of ages of between 2.5 - 5 Myr. Such a conclusion is supported 
by inspection of the properties of individual stars in both associ- 
ation and cluster, although interpretation is complicated by pos- 
sible contamination of the latter by the former if, as seems likely, 
Sco OB 1 hosts a younger population than is found in NGC 623 1 . 

A large number of Main Sequence (MS) stars are found 
within NGC 6231, with the 08 V companion in the binary 
HD 152234 apparently being the earliest and defining the MS 
turnoff; consistent with an age of 5 Myr. While an 06 V compan- 
ion to the WC7 star WR79 has been reported (Hill et al. 2000) 
the 126 day binary period of HD 152234 will have ensured that 
neither component in this system will have interacted, while the 
8.89 day period of WR79 indicates that significant mass trans- 
fer to the secondary may have occured (e.g. Petrovic et al. 2005). 
Encouragingly, the 09.71a primary of HD 152234 also lies on the 
5 Myr isochrone for rapid initial rotation, with the current vsim 
also being consistent with this placement (Fraser et al. 120101 . 

Building on this approach, we find that HD 152219 (09.5 
III + Bl-2 III-V) and HD 152134 (B0.5 la) lie upon the 5 Myr 
evolutionary track for non-rotating stars, with - sequentially - 
HD 326331 (08 III((f))), HD 152247 (09 III), HD 152249 (09 
Ib((f))) and HD 152234 (09.7 la + 08 V) following the rotating 
track. In both cases the systematic progression to later spectral 
types with increasing luminosity class suggest these are real evo- 
lutionary sequences (cf. Cyg OB2; Negueruela et al. 2008). 

However, there are a number of stars within NGC 623 1 and 
Sco OB 1 that appear incompatable with a ~ 5 Myr population: 
the 08 Iafpe/WN9ha stars HD 151804 and 152048 and the 05.5- 
7.5 III stars HD 151515, 152723, 152233 and 15 2248. Regarding 
the former pair, Bohannan & Crowther d 19991 1 suggest a close 
physical kinship between these objects, apparently precluding an 
evolution through the cooler late-O/early-BSG phase present in 
NGC 6231. Of the mid-O giants, HD 151515 and 152723 appear 
somewhat subluminous for their spectral type but the properties 
of the 05.5 m(f) + 07.5 and 07.5 111(f) + Q7 1 11(f) binaries 
HD 152233 and 152248 (Sana et al. I20TJ81 [200TT) clearly sup- 
port the assertion that a younger population is present; although 
the dymanical masses for both components of the latter system 
appear somewhat lower than expected (29. 6M© and 29.9M re- 
spectively). 

We therefore conclude that comparison of NGC 6231 to 
Sco OBI implies that the respective stellar populations are 
non coeval. However, with the exception of HD 152233 and 
HD 152248, the massive stellar population of NGC623 1 appears 
consistent with a single burst of star formation ~5 Myr ago. 
Indeed, no early-mid O MS stars consistent with a younger pop- 
ulation appears present within NGC 6231, while the population 
of late-O/early-B MS stars that are present are systematically 
displaced redwards from the <3 Myr isochrones. Therefore, if 
they are members of NGC 6231, we cannot exclude the possi- 
bility that, with spectral types <08, HD 152233 and HD 152248 
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are in fact bona fide blue stragglers; having evolved via a differ- 
ent pathway from the majority of stars and hence that NGC 623 1 
is truly coeval. In this respect it would therefore closely resem- 
ble Cyg OB2, for which Negueruela et al. (2008 ) arrived at a 
comparable conclusion. 

Appendix D: Spectral fits to ^ Sco and HD1 90603 

We present additional model fits to f Sco based on the FEROS 
(3700-8850A) optical data (not presented in the paper version 
due to reasons of space). 
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Fig. D.l. Comparison of the synthetic spectrum of ^' Sco (red dashed line) to observational data for various transitions between 
4090-22 100A. Optical data correspond to the FEROS data kindly provided by Otmar Stahl. Note the lack of long-termvariability 
between these data and those presented in Fig. 7. 
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Fig.D.2. Same as Fig ED Cont. 
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